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ABSTRACT 

Introductory lectures on heavy quarks and heavy quark effective field 
theory. Applications to inclusive semileptonic decays and to interactions 
with light mesons are covered in detail. 

1. Introduction 

The theory of heavy hadrons has received much attention over the last few 
years. There are two good reasons for this. On the one hand experiments probing 
the properties of heavy hadrons have come of age. Not a year goes by without 
some remarkable new measurement, and 1994 is no exception[l]. Moreover, SLAC 
is moving furiously towards completion of a B factory, which will hopefully operate 
well before the end of the decade. On the other hand, there have been many 
theoretical advances that allow clearer interpretation of the experimental results. 
These lectures are an introduction to these theoretical developments. 

The lectures are not intended to be encyclopedic in any one subject. I have 
decided to try to convey the principal ideas as clearly as I can, and to give some 
sample applications here and there. Occasionally I may describe "the state of the 
art" in a given field, without necessarily entering into details. I hope to have 
included enough references that the reader may follow up on any of these subjects 
if so inclined. 

Exercises are scattered throughout the lectures. In the age of electronic typog- 
raphy it was easy enough to display the exercises in smaller print and separate them 
with clearly visible horizontal lines. I hope that the material given in the lectures 
is sufficient for obtaining the solution to the problems. 

Chapter 2 contains a brief description of what Effective Field Theories are, at 
least in the very specific context of weak interactions. The presentation is somewhat 
telegraphic. I expect the student to know something about effective lagrangians. 



The intention is to show you how I think about the subject so that the presentation 
of the effective lagrangian of heavy quarks, in Chapter 3, goes down more easily. 

Chapter 3 describes the Heavy Quark Effective Field Theory (HQET) and its 
symmetries, to leading order in the large mass. Some examples and applications 
are given. The corrections of order 1/M are described in Chapter 4. 

The last two Chapters concentrate on applications of the HQET that have 
received a lot of attention over the last year. This is where these lectures deviate 
substantially from my Mexican lectures. Chapter 5 gives the proof that the rate for 
inclusive semileptonic B decay is given by the parton decay rate, while Chapter 6 
introduces an effective lagrangian of heavy mesons and pion, interacting in a chirally 
invariant way, and respecting heavy quark symmetries. 



2. Preliminaries 

2.1. Conventions and Notation 

The metric is (H ). Gamma matrices satisfy {7^, 7^} = 2g fJil/ , 7 is hermi- 

tian, 7* antihermitian, 75 = i7°7 1 7 2 7 3 and <7 Mi , = -§[7^, 7"]. In the Dirac convention 
7 = diag(l, 1, -1, -1). e 0123 = +1. 

Left and right handed fields are denoted by subscripts, 

i>L = |(i - 75)^ i>R = 5(1 + 75)^ 

States have the relativistic normalization, 

(p\p') =2E5(p-p') (2.1) 

unless otherwise noted. 

The charged current interaction in the standard model is given by the lagrangian 

C int = ^=W^(u c £)7m(1-7s)^ ( s ] + h.c. , (2.2) 
where V is the 3x3 unitary CKM matrix. 



2.2. Effective Lagrangians 

The study of decays of heavy mesons involves two very different types of physics. 
On the one hand there is the underlying interactions that are responsible for the 
decay. These may be the ordinary weak interactions of the standard model or 
some new, presumably weaker, interactions, e.g., extended technicolor or exchange 
of scalar quarks in supersymmetric theories. On the other hand there are strong 
interactions which modulate the rates of the decays. Ultimately we would like to 
uncover the former, which requires some level of understanding of the latter. 

Since the two physical aspects are conceptually different, it seems natural to 
separate them in our computations. This is where effective lagrangians come in 
handy. An effective lagrangian for meson decays will only involve the dynamical 
degrees of freedom that are relevant. For example, light quarks (it, d, s) and gluons 
for K meson decays. The interaction responsible for the decay, a W exchange, is 
represented as a AS = 1 four-quark operator. The fact that the TU-boson is no 
longer in the theory does not impair our ability to predict the decay rate, up to an 
accuracy or order m^/M^. 

Is this a major step backwards? After all, you may argue, this is going back to 
the Fermi theory of weak interactions. Effective lagrangians prove useful because, 

(i) The computation of long distance physics is decoupled from the short distance 
physics. Thus, one can make a catalog of matrix elements of operators between 
meson states. This could then be used to compute the effects of any fundamental 
theory, after reducing its lagrangian to an effective one. 

(ii) They provide a method for the computation of amplitudes when disparate scales 
are present. In iiT-meson decays the ratio m^/Mjy is a small number. Loga- 
rithms of this ratio invalidate even the perturbative computation of the decay. 
Effective lagrangians provide a method for resummation of these large logs. 

(iii) One can characterize, d la Fermi, all possible interactions in terms of operators. 
In other words, one can make model independent analysis of the possible effects 
of new physics. 

(iv) It's the right way to think about the low energy effects of very heavy particles. 

2.3. Formulating Effective Lagrangians 

The precise meaning of effective lagrangians is best formulated in terms of 
relations between Green functions. Take again the example of weak interactions 
at low energies, that is, when all the momenta involved are much smaller than the 



VF-boson mass. Everyone knows that we can account for the effects of the VF-boson 
by adding to the Lagrangian terms of the form 

A£ eff = -t^kO , (2.3) 

IVI W 

where O is a 4-fermion operator and k contains mixing angles and factors of the 
weak coupling constant. This is simply the statement that a Green function G of 
the original theory (the standard model including QCD) can be approximated by a 
Green function Go of the effective theory (a gauge theory of QCD and electromag- 
netism) with an insertion of the effective Lagrangian: 

G=^-kGo + ... . (2.4) 

The ellipses stand for terms suppressed by additional powers of [M w )~ 2 . This 
equation replaces the task of computing the more complicated left side, which de- 
pends on Mw, by the computation in the effective theory which is independent of 
Mw, and indeed, completely free of the VF-boson dynamical degrees of freedom. 
On the right hand side, the factor of l/M w gives the dependence on the W-boson 
mass. 

The full theory has logarithmic dependence on Mw which has not been made 
explicit. Eq. (2.4) is not quite correct. The correct version is[2] 

G=^KC(M w /^g s )G + ... • (2.5) 

The function C is, in this case, also known as the 'short distance QCD effect' first 
calculated for the process s — > uud by Altarelli and Parisi[3], and Gaillard and 
Lee [4]. 

Summing up, an effective theory (of either the 'normal' or the HQ type) is a 
method for extracting explicitly the leading large mass dependence of amplitudes. 
Moreover, the rules of computation of the effective theory are completely indepen- 
dent of the large mass. 

2.4- Computing Effective Lagrangians 

At tree level the computation of effective lagrangians is straightforward since 
the short distance QCD effects can be neglected. This is best explained through 
an example. Consider transitions with AC = —AU = AS = —AD = 1. In the 
full theory these take place by the exchange of a VF-boson coupling to the charged 



currents (clI^sl) and {(Il^^ul)- Take the Green function for, say, c — > swci and 
expand in powers of momentum over the VF-boson mass. This amount to expanding 
the VF-boson propagator; in 'tHooft-Feynman gauge, 
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pi - Ml 

The effective hamiltonian for the decay is then 

ft eff = ^v: d V cs {c Ll »s L ){d Ll »UL) + ■■■ (2.6) 

where we have introduced Fermi's constant Gp = \f2g\j%Ml. The ellipsis stand 
for operators of higher dimension which come from the expansion of the propagator, 
replacing p — > id, as in 

(c L rs L )d 2 {d L Yu L ) 

Exercise 2.1 Show that the effective hamiltonian for AB = 1, AC = AU = 0, is given at 
tree level by 

HeS = ^fJ2 J2 KbV qq '(bL a 7 fl qL a )(qL^q' L(} ) (2.7) 

q = u,c q> =d,s 

where a and (3 are color indices. 

Beyond tree level the simple effective hamiltonian of Eq. (2.6) is replaced by 
a sum over operators of dimension six and with the same quantum numbers. Ne- 
glecting the masses of the light quarks there is only one more operator: 

H eS = ^=-V* d V cs CiO t + ■■■ (2.8) 

where 

1 = (c Ll ^s L )(d Ll ^u L ) 

2 = {c L ^u L ){d L ^s L ) 

and the ellipsis stand again for higher dimension operators suppressed by additional 
powers of Gf- The coefficients q encode all the information about the dependence 
on the W-mass (beyond the trivial factor of Gf)- We already know that at tree 
level ci = 1 and = 0. They may be computed as an expansion in a s . It is best, 



however to reorganize the expansion to account for the large ratio of scales m c /M w 
by assuming a s C 1 and a s \og(m c /Mw) ~ 1, thus 
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Here 6o = 25/3 stands for the coefficient of the first term in the perturbative expan- 
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sion of the beta function in QCD, (3(g) = — bo + . . ., in the case of four flavors 
of quarks. This is the so-called "leading-log" approximation to the coefficients. We 
have displayed order of the "sub-leading-log" or "next-to-leading-log" corrections. 

Exercise 2.2 Write down the complete list of operators that contributes to the effective 
hamiltonian for AB = 1, AC = AU = 0. Neglect the masses of u, d, s and c-quarks, but 
not of the 6-quark. (Hint: What is the symmetry group of QCD in the massless limit? 
How do the operators in Eq. (2.7) transform under this symmetry?) 

The computation of the coefficients is itself quite simple but it detracts from 
our main focus. A rather similar computation is presented below for the Heavy 
Quark Effective Theory; see Section 3.9. For a beautiful exposition of the method, 
see [2]. 

Exercise 2.3 When you understand the procedure described in Section 3.9 return here 
and check Eqs. (2.9). What modifications are necesary to account for a fifth quark with 
mass nib intermediate between charm and the weak scale? (Assume m c <C m& <C Mw)- 



Exercise 2.4 Show that the effective hamiltonian for AB = -AD = 1, AC = AU = of 
the previous exercises can be written as the sum of exactly two terms, 

n w = ^(CcO c + CtO t ) (2.io) 

V2 

where £ g = V* b V q d, and O q are linear combinations of composite operators that may 
depend on Mw but not on CKM angles. Use phenomenological information to show that 
in the effective hamiltonian for AB = —AS = 1, AC = AU = 0, depends only on the 
single mixing angle V^Vcs- 



3. Heavy Quark Effective Field Theory 

3.1. Intuitive Introduction 

The central idea of the HQET is so simple, it can be described without reference 
to a single equation. And it should prove useful to refer back to the simple intuitive 
notion, to be presented below, wherever the formalism and corresponding equations 
become abstruse. 

The HQET is useful when dealing with hadrons composed of one heavy quark 
and any number of light quarks. More precisely, the quantum numbers of the 
hadrons are unrestricted as far as isospin and strangeness, but are ±1 for either 
B- or C- number. In what follows we shall (imprecisely) refer to these as 'heavy 
hadrons'. 

The successes of the constituent quark model is indicative of the fact that, inside 
hadrons, strongly bound quarks exchange momentum of magnitude a few hundred 
MeV. We can think of the typical amount A by which the quarks are off-shell in 
the nucleon as m p /3 330MeV. In a heavy hadron the same intuition can be 
imported, and again the light quark(s) is(are) very far off-shell, by an amount of 
order A. But, if the mass Mq of the heavy quark Q is large, Mq 3> A, then, in 
fact, this quark is almost on-shell. Moreover, interactions with the light quark(s) 
typically change the momentum of Q by A, but change the velocity of Q by a 
negligible amount, of the order of A/Mq <C 1. It therefore makes sense to think of 
Q as moving with constant velocity, and this velocity is, of course, the velocity of 
the heavy hadron. 

In the rest frame of the heavy hadron, the heavy quark is practically at rest. 
The heavy quark effectively acts as a static source of gluons. It is characterized 
by its flavor and co\or-SU(3) quantum numbers, but not by its mass. In fact, 
since spin-flip interactions with Q are of the type of magnetic moment transitions, 
and these involve an explicit factor of g s /MQ, where g s is the strong interactions 
coupling constant, the spin quantum number itself decouples in the large Mq case. 
Therefore, the properties of heavy hadrons are independent of the spin and mass of 
the heavy source of color. 

The HQET is nothing more than a method for giving these observations a 
formal basis. It is useful because it gives a procedure for making explicit calcula- 
tions. But more importantly, it turns the statement 'Mq is large' into a systematic 
perturbative expansion in powers of A/Mq. Each order in this expansion involves 
QCD to all orders in the strong coupling, g s . Also, the statement of mass and spin 



independence of properties of heavy hadrons appears in the HQET as approximate 
internal symmetries of the Lagrangian. 

Before closing this Section, we point out that these statements apply just as 
well to a very familiar and quite different system: the atom. The role of the heavy 
quark is played by the nucleus, and that of the light degrees of freedom by the 
electrons (and the electromagnetic field) 1 . That different isotopes have the same 
chemical properties simply reflects the nuclear mass independence of the atomic 
wave- function. Atoms with nuclear spin s are 2s + 1 degenerate; this degeneracy is 
broken when the finite nuclear mass is accounted for, and the resulting hyperfine 
splitting is small because the nuclear mass is so much larger than the binding 
energy (playing the role of A). It is not surprising that, using Mq independence, 
the properties of B and D mesons are related, and using spin independence, those 
of B and B* mesons are related, too. 

3.2. The Effective Lagrangian and its Feynman Rules 

We shall focus our attention on the calculation of Green functions in QCD, 
with a heavy quark line, its external momentum almost on-shell. The external 
momentum of gluons or light quarks can be far off-shell, but not much larger than 
the hadronic scale A. This region of momentum space is interesting because physical 
quantities — S'-matrix elements — live there. And, as stated in the introduction, 
we expect to see approximate symmetries of Green functions in that region which 
are not symmetries away from it. That is, these are approximate symmetries of a 
sector of the S'-matrix, but not of the lagrangian. 

The effective Lagrangian C e s is constructed so that it will reproduce these Green 
functions, to leading order in A/Mq. It is given, for a heavy quark of velocity 
(v 2 = 1), by[5], 

4ff = Qviv-DQ V , (3.1) 

1 An obvious distinction between the atomic and hadronic systems is that in the latter 
the configuration of the light degrees of freedom is non-computable, due to the difficulties 
afforded by the non-perturbative nature of strong interactions. The methods that we are 
describing circumvent the need for a detailed knowledge of the configuration of light degrees 
of freedom. The price paid is that the range of predictions is restricted. To emphasize 
the non-computable aspect of the configuration of light degrees of freedom, Nathan Isgur 
informally referred to it as "brown muck", and the term has somewhat made it into the 
literature. 



where the covariant derivative is 

Dp = <9 M + ig s A®T a , (3.2) 

and the heavy quark field Q v is a Dirac spinor that satisfies the constraint 

1 + A 

Qv = Qv ■ (3-3) 

In addition, it is understood that the usual Lagrangian £n g ht for gluons and light 
quarks is added to C^. 

We can see how this arises at tree level, as follows [6]. Consider first the tree 
level 2-point function for the heavy quark 

We are interested in momentum representing a quark of velocity slightly off-shell: 

p M = Mqv^ + fc M . (3.5) 

Here, 'slightly off-shell' means is of order A, and independent of Mq. Sub- 
stituting in Eq. (3.4), and expanding in powers of A/Mq, we obtain, to leading 
order, 



v-k V M, 



We recognize the projection operator of Eq. (3.3), and the propagator of the lagran- 
gian in (3.1). 

Similarly, the 3-point function (a heavy quark and a gluon) is given by 

G ? '"•<» = ^-^nj^-^M, (3.7) 

where A vlx (q) is the gluon propagator. Expanding as above, we have 

G ; w, ^»=( 1 i 1 )s(-» r "^ A ^'» +o (4)' (3 - 8) 



where we have used 



Again, this corresponds to the vertex obtained from the effective Lagrangian in 
Eq. (3.1). 



Exercise 3.1 Extend these results to arbitrary tree-level Green functions (but only those 
with one heavy quark and all other (light) particles carrying momentum of order A). 

The effective Lagrangian in (3.1) is appropriate for the description of a heavy 
quark, and indeed a heavy hadron, of velocity v^. It does, however, break Lorentz 
covariance. This is not a surprise, since we have expanded the Green functions 
about one particular velocity: in boosted frames, the expansion in powers of A/Mq 
becomes invalid, since the boosted momentum k^ can become arbitrarily large. 
Lorentz covariance is recovered, however, if we boost the velocity 

-> Apvv,, (3.10) 

along with everything else. It will prove useful to keep this simple observation in 
mind 2 . 

Just as in the case of an effective theory for weak interactions, when one goes 
beyond tree level one must be careful to make explicit any anomalous mass depen- 
dence. When the W is integrated out the correction was encrypted in some function 
C(M\v) in Eq. (2.5). The situation is entirely analogous in the HQET. We have 
introduced an effective Lagrangian such that Green functions G v (k;q) calcu- 
lated from it agree, at tree level, with corresponding Green functions G(p; q), in the 
QCD to leading order in the large mass 

G{p; q) = G v (k; q) + O (A/Mq) (tree level) . (3.11) 

Here, A stands for any component of /c M or of the q's, or for a light quark mass, and 
p = Mqv + k. We will come back to the study of the HQET beyond tree level and 
will make explicit the anomalous mass dependence in Section 3.8. 

2 In an alternative method, championed by Georgi[7], the effective Lagrangian C e s 
consists of a sum over the different velocity Lagrangians, 

4ff, of Eq. (3.1). Lorentz 

invariance is recovered at the price of "integrating in" the heavy degrees of freedom. This 
does not lead to overcounting of states, because the sectors of different velocity do not 
couple to each other, a fact that Georgi refers to as a "velocity superselection rule" . See 
also [8]. 



3.3. Symmetries 

Flavor - SU(N) The Lagrangian for N species of heavy quarks, all with velocity 

v, is 

N 

4? = £<5?> iv - D • ( 3 - 12 ) 

i=i 

This Lagrangian has a U(N) symmetry [9, 10, 11]. The subgroup U(1) N corresponds 
to flavor conservation of the strong interactions, and was a good symmetry in the 
original theory. The novelty in the HQET is then the nonabelian nature of the 
symmetry group. This leads to relations between properties of heavy hadrons with 
different quantum numbers. Please note that these will be relations between hadrons 
of a given velocity, even if of different momentum (since typically Mq. ^ Mq. for 
i 7^ j). Including the b and c quarks in the HQET, so that N = 2, we see that the 
B and D mesons form a doublet under Qavor-SU(2). 

This flavor-5 , L ? "(2) is an approximate symmetry of QCD. It is a good symmetry 
to the extent that 

m c > A and m b > A . (3.13) 

These conditions can be met even if m^ — m c 3> A. This is in contrast to isospin 
symmetry, which holds because rrid — m u <^ A. 

In the atomic physics analogy of Section 3.1, this symmetry implies the equality 
of chemical properties of different isotopes of an element. 

Spin - SU(2) The HQET Lagrangian involves only two components of the 
spinor Q v . Recall that 

1 ^ JL jQv = 0. (3.14) 

The two surviving components enter the Lagrangian diagonally, i.e., there are no 
Dirac matrices in 

4ff = Qv ivD Q v . (3.15) 

Therefore, there is an SU(2) symmetry of this Lagrangian which rotates the two 
components of Q v among themselves [9, 12-13]. 

Please note that this "spin" -symmetry is actually an internal symmetry. That 
is, for the symmetry to hold no transformation on the coordinates is needed, when 
a rotation among components of Q v is made. On the other hand, to recover 
Lorentz covariance, one does the usual transformation on the light-sector, including 
a Lorentz transformation of coordinates and in addition a Lorentz transformation 
on the velocity ly. A spin-5'L r (2) transformation can be added to this procedure, 
to mimic the original action of Lorentz transformations. 



Exercise 3.2 To make it plain that this symmetry has nothing to do with "spin" in the 
usual sense, consider the large mass limit for a vector particle [14]. Use the massive vector 
propagator 

_ z ^-w^ 2 (316) 

p 2 — m A 

to obtain the Lagrangian for the HVET (Heavy Vector Effective Theory) 

4£ = 4„ iv-DA vlx , (3.17) 

with the constraint 

{v^v v - g^ v )A vv = A v>1 . (3.18) 

What is the dimension of this effective vector field? Why? Show that the effective La- 
grangian is invariant under an SU (3) group of transformations, rotating the three com- 
ponents of the vector field among themselves. Note that the "spin" symmetry is not 
associated with SU(2) in this case. 

The symmetry of the theory is larger than the product of the flavor and spin 
symmetries. If there are Ns, Np, and Ny species of heavy scalars, fermions, and 
vectors, respectively, all transforming the same way under coior-5'L r (3), the sym- 
metry of the effective theory is SU(N S + 2N F + 3N V ). 

Exercise 3.3 What is the symmetry group for a theory with Ns, Np, and Nv species 
of heavy scalars, fermions, and vectors, respectively, all transforming the same way under 
color-S'C/(3), in D space-time dimensions? (If you can't handle arbitrary D, try D = 2 
and D = 3). 



3.4- Spectrum 

The internal symmetries of the effective Lagrangian are explicitly realized as 
degeneracies in the spectrum and as relations between transition amplitudes. In 
this Section we will consider the spectrum of the theory [15]. 

Keep in mind that momenta, and therefore energies and masses, are measured 
in the HQET relative to Mqv^. Therefore, when we state that in the HQET the B 
and D mesons are degenerate, the implication is that the physical mesons differ in 
their masses by — m c . 

For now let us specialize to the rest frame v = (1,0). The total angular mo- 
mentum operator J, i.e., the generator of rotations, can be written as 

J = L + S, (3.19) 

where L is the angular momentum operator of the light degrees of freedom, and S, 
the angular momentum operator for the heavy quark, agrees with the generator of 



spin-SU(2). Since J and S are separately conserved, L is also separately conserved. 
Therefore, the states of the theory can be labeled by their L and S quantum numbers 
(I, mf, s,m s ). Of course, s = 1/2, so m s is 1/2 or —1/2 only. 

The simplest state has / = and, therefore, J = 1/2. We will refer to it as the 
Aq, by analogy with the nonrelativistic potential constituent quark model of the 
A-baryon, where the strange quark combines with a / = 0, I = 0, combination of 
the two light quarks. 

Next is the state with / = 1/2. It leads to J = and J = 1. We deduce that 
there is a meson and a vector meson that are degenerate. For the 6-quark, the B 
and B* fit the bill. They are the lowest lying B = — 1 states. The lowest lying 
C = 1 states are the D and D* mesons. These again can very well be assigned 
to our J = and J = 1 multiplet. The difference M D * - M D = 145 MeV is 
reasonably smaller than the splitting between the D* and the next state, the D\, 
with M Dl - M D * = 410 MeV. 

The splittings of B and B* and of D and D* result from spin-SU(2) symmetry 
breaking effects. These must be corrections of order A/Mq to the HQET predic- 
tions. Therefore, one must have Mb* — Mb = A 2 /mj and analogously for the D-D* 
pair. Therefore 

Mb* - M B = 
M D * - M D m b 

Approximating m c and mj, by Mb, and Mb, respectively, we get ~ 1/3 on the 
right side, in remarkable agreement with the left side. Although these results also 
follow from potential models of constituent quarks, it is important that they can be 
derived in this generality, and this simply. 

The states with / = 3/2 have J = 1 and 2. The D\ and D\, with Md* — 
= 40 MeV, are remarkably closely spaced (and of course, have the appropriate 
quantum numbers to form a spin multiplet). 

Exercise 3.4 A more complete classification of the spectrum would include parity. What 
modifications are needed, if any, to include parity as a quantum number? You should find 
that there are two possible £ = 1/2 multiplets. One corresponds to the D and D* (or 
B and B*) mesons. It is usually argued that the multiplet with opposite parity would 
contain very broad resonances that could not be identified as stable states. Why? 

While in the infinite mass limit states \l,mi;s,m 8 ) have sharp L 2 , L Z: S 2 and S z , 
these are not good quantum numbers for physical states. Regardless of how small 
spin-symmetry breaking effects may be, they force states into linear combinations 
of sharp J 2 , J 2 , L 2 and S 2 , \ J,mj;l, s). SU (2)-spin transformations connect states 
of J = / + 1/2 with those of J = / - 1/2. Now 

\J,mj;l,s) = ^\l,mi;s,m 8 )CJ™ J 8ms 



where C/^ J sms = C(lmi; sm s \Jmj) are Clebsch-Gordan coefficients. The decom- 
position is useful because we know how the states on the right transform under 
spin-5 , L r (2). The inverse expression, 

|Z, m,; s, m s ) = ^ | J, m j; /, s) (Cf™? ams )* 

gives the linear combinations of physical states with definite spin- ST/ (2) numbers. 

For example, for the B and B* multiplet, the mi = 1/2 and mi = —1/2 states 
that form spinet/" (2) doublets are, respectively 

/ x / iT(0)-B x 

^1/2= (^b^o^bJ and ^-1/2= ^ b*{-) ) ' 

Rotations mix components among these doublets. We can combine them into a 
matrix \1/ Qa = (ip a )a- If T>^(R) stands for a 21 + 1 dimensional representation of 
the rotation i?, then the action of spin-SU(2) alone is ^ — > "D^ 1 / 2 ^/?)^/, while a 
rotation is ^ -> p( 1 /2)( jR )^p(i/2)( jR )t. 

This is easily generalized. For arbitrary / there are 2/ + 1 doublets of spin- 
SU{2), (f) a , a = —/,..., /. They can be assembled into a 2 x (2/ + 1) matrix <E> aa = 
(0a)cn which transforms as $ — > T'( 1 / 2 )(i?)$r'^)(i?) t under rotations. The linear 
combination of physical states in $ aa can be written as a sum of at most two terms: 

$ =y(+) A + y(-) A 

where (x«o ) is the state with J = I + 1/2 (J = / — 1/2), A = mj = a + a, 

weighted by the corresponding Clebsch-Gordan coefficient, C^y^a • 

If Qo is the two component heavy quark field for v = (1,0) then, according to 
the Wigner-Eckart theorem the matrix elements of QqTQq, for any Pauli matrix T, 
between B and B* states are all given in terms of a single reduced matrix element, 
X, times appropriate symmetry factors constructed out of Clebsch-Gordan coeffi- 
cients. Let the state = J2 aa ^aa|«a). Then the matrix element (iP'\QqTQo\iP) 
must be linear in ^ aa -, i $>'^ b and r*, and proportional to the invariant tensor 
{l3b\Q^ Qos\aa) . There is only one invariant one can construct out of these tensors. 
This can be summarized as follows, 

<<//|QjrQ |V'>=xM'rv|>. (3.22) 

By this we mean that if the state ip is, say, a S-meson, then the corresponding 
matrix \1/ on the right hand side is obtained from the matrix \l/ aa of Eq. (3.21) by 
setting B* = and 5 = 1, and analogously for the other possible choices of the 
states ip and ift'. In the next Section we generalize this result to the case were the 
states if) and tfj' may have different velocities. 



3. 5. Covariant Representation of States 

In the Chapters that follow we will be interested in extracting the consequences 
of the spin and flavor symmetries of the HQET to a variety of processes. These 
processes may involve transitions between heavy hadrons of different velocities. It 
is convenient to develop a formalism that automatically extracts the information 
encoded in the symmetries [16]. I follow the simple presentation of Ref. [17]. 

A prototypical example of an application is the computation of relations be- 
tween form factors in semileptonic B to D and D* decays. There one needs to study 
the matrix elements 



where uq is a spinor representing the heavy quark, t/juq = uq, and v q is an antispinor 
representing the light stuff with / = 1/2, satisfying v q f = v q . The product in (3.24) 
is a superposition of states with J = and 1. To identify the pseudoscalar meson P 
and the vector meson V(e) with polarization e, ev = 0, we must form appropriate 
linear combinations of the spin up and down spinors. This is most easily done in 
the rest frame v = (1,0); the result will be generalized to arbitrary v by boosting. 
In the Dirac representation the spin operator is S = 7 5 7°7/2 so that the spinor 
basis Ua = o"ia an d u«' = #2a corresponds to spin up and spin down, and the 
antispinor basis Va = ~ ^3a and = — &4 Q corresponds to spin down and spin 
up. With S(uv) = (Su)v + u(Sv) it is easy to check that the combination 



(D(v)\cvTb v >\B(v')) and (D*(v),e\c v Tb v >\B(v')) . 



(3.23) 



We would like to represent these / = 1/2 mesons as the product 



UQV q , 



(3.24) 




(3.25) 



has zero spin, while 




(3.26) 



with = (0, l,±z,0) and e^ ) = (0,0,0,1), have total spin 1, with third com- 
ponent 1, and —1, respectively. Thus, for arbitrary velocity v one obtains the 
representation for pseudoscalar and vector mesons: 

M W= (i±i) 7 s M' M =(±±t)j. (3.27) 

By construction, the spin symmetry acts on this representation only on the first 
index of the matrices M(v) and M*(v, e). 

The power of this machinery can now be displayed. Consider the matrix ele- 
ments (3.23). Using the above representation of states and noting that the result 
should transform under the spin symmetry just as the matrix T, we have 

(D(v)\c v Tb v >\B(v')) = -£(vv')TrD(v)TB(v') (3.28a) 
(D*(v)e\c v Tb v ,\B(v')) = -£(W)Tr D*(v, e)TB{v') , (3.286) 

where X = 7°X^7°. The common factor — £(w') plays the role of the reduced 
matrix element in the Wigner-Eckart theorem. We will explore the consequences 
of Eqs. (3.28) in depth in Section 3.7. 

Exercise 3.5 An even simpler case is that of the I = multiplet. In this case the states 
must transform as a spinor. How would you represent matrix elements of these states? 
What about those of the I = 1 or I = 3/2 multiplets? This formalism can be extended[17] 
to deal with multiplets of arbitrary I. 

3.6. Meson Decay Constants 

The pseudoscalar decay constant is one of the first physical quantities studied 
in the context of HQET's. For a heavy-light pseudoscalar meson X of mass Mx, 
the decay constant fx, we will see, scales like 1 / V Mx ■ This was known before the 
formal development of HQET's, although the arguments relied on models of strong 
interactions. The HQET will give us a systematic way of obtaining this result. 
Moreover, it will give us the means of studying corrections to this prediction. 

The decay constant fx is defined through 

(0\A^0)\X(p)) = f x p^ , (3.29) 

where = qj^jsQ is the heavy-light axial current, and the meson has the standard 
relativistic normalization 



(X{p')\X(p)) = 2Ed^(p-p') 



(3.30) 



Thus, the states have mass-dimension — 1. Analogous definitions can be made for 
other mesons. For example, for the vector meson X* (the 1 = 1/2 partner of X), 
has 

(0|y M (0)|X*(p,e)> = /x*e M . (3.31) 

Note that the mass-dimensions of fx and fx* are 1 and 2, respectively. 

Consider the decay constant of the meson state in the HQET. The effective 
pseudoscalar decay constant fx is defined by 

(0\A^(0)\X(v)) = f x v„ (3.32) 

The state in the HQET, \X), is normalized a la Bjorken and Drell,[18] to 2E/M X 
rather than to 2E: 

(X(v')\X{v)) = 2v°5W(v - v') . (3.33) 

Actually, defining states in the HQET requires some care, but I will just assume 
it all works and merely refer the interested reader to the literature [8]. Obviously, 
since the normalization of states and the dynamics are Mq independent, so is fx- 
To relate fx to the physical fx simply multiply Eq. (3.32) by \fWx, to restore the 
normalization of states of Eq. (3.30), and write = p^/Mx- Thus we arrive at 



fx = fx/VMx, (3.34) 
A useful way of quoting the result is, for the physical case of B and D mesons, 



As a simple application of the spin symmetry, consider the pseudoscalar decay 
constant fx*- Using the 4x4 notation of Section 3.5, the matrix element in 
Eq. (3.31) that defines the pseudoscalar constant is proportional to 

Tr ( 7 ^7 5 M(t;)) = Tr (V 75 (^) 7s) = (3.36) 
The matrix element 

(0\V»(0)\X*(v)e) = f x *e» (3.37) 

is proportional to 




(3.38) 



with the same constant of proportionality. Therefore 

fx* = -fx (3.39) 
The sign is unimportant, since it can be absorbed into a phase redefinition of either 
state. It is the magnitude that matters. Multiplying by \[M~x* ~ \[Mx to restore 
to the standard normalization, we have 

f x * = -f x M x (3.40) 

The predictions Eq. (3.35) and Eq. (3.40) have not been tested experimentally. 
The difficulty is the small expected branching fraction for the decays X — > \xv or 
X* — > for X = B and D. Alternatively, the decay constants fx and fx* can be 
measured in Monte Carlo simulations of lattice QCD. There are indications from 
such simulations that the 1/Mq corrections to the relation (3.35) are large [19]. 

3.7. Semileptonic decays 

The semileptonic decays of a B- meson to D- or D* -mesons offer the most direct 
means of extracting the mixing angle \V c b\. In order to extract this angle from 
experiment, theory must provide the form factors for the B — > D and B — > D* 
transitions. Several means of estimating these form factors can be found in the 
literature. A popular method consists of estimating the form factor at one value of 
the momentum transfer q 2 = q 2 , and then introducing the functional dependence 
on q 2 in some arbitrary, hopefully reasonable, way. In the pre-HQET days it was 
customary to estimate the form factor at Op from some model of strong interactions, 
like the non-relativistic constituent quark model. 

The HQET gives the form factor at the maximum momentum transfer, q 2 = 
(/max = (Mb — Md) 2 — the point at which the resulting D or D* does not recoil in 
the restframe of the decaying S-meson. While the functional dependence on q 2 is 
a non-perturbative problem, it is already progress to have a prediction of the form 
factor at one point. Moreover, the HQET gives relations between the form factors. 
One may study these relations experimentally to test the accuracy of the HQET 
predictions. 

The standard definition of form factors in semileptonic S-meson decays is 
(D(p')\V^\B(p)) = f + (q 2 )(p + p')^ + f_(q 2 )(p-p') fl (3.41a) 
{D*(p')e\A^\B(p)) = /(g 2 ) e ; + a + (q 2 y-p(p + p% + a-(q 2 )e*-p(p - p% (3.416) 

(D*(p')e\V^\B(p))=zg(q 2 )e^ Xa e^(p + p') x (p-pT (3.41c) 
Here, the states have the standard normalization, Eq. (3.30), and q 2 = (p — p') 2 . 
The contribution to the decay rates from the form factors /_ and a_ are suppressed 
by m 2 /M B , where me is the mass of the charged lepton, and therefore they are often 
neglected. 



Exercise 3.6 Compute the differential decay rates, dT/dx dy, where x = q 2 /Mg and 
y = v ■ q/Ms, for B — > Dei? and B — > D*ev, in terms of these form factors. 

In the effective theory, we would like to compute the matrix elements of the 
effective currents and between states of the / = | multiplet. We can take 
advantage of the flavor and spin symmetries to write these matrix elements in terms 
of generalized Clebsch-Gordan coefficients and reduced matrix elements, i.e., we use 
the Wigner-Eckart theorem. We have already introduced the relevant machinery in 
Section 3.5. The matrix elements of the operator G = c v >Tb v between B and D or 
D* states, are given by (c.f., Eqs. (3.28)) 

(D(v')\G\B(v)) = -C(v-v')TtD(v')TB(v) (3.42a) 
(D*(v')e\G\B(v)) = -£(vv')TrD*(v',e)rB(v). (3.426) 

Before expanding Eqs. (3.42), we note that the flavor symmetry implies that 
the S-current form factor between S-meson states is given by the same reduced 
matrix element: 

(B(y')\b v >rb v \B(y)) = -£(W)Tr B(v')TB(v) (3.43) 

Using r = 7 , and recalling that B- number is conserved, one finds that £ is fixed 
at d' = v. With the normalization of states appropriate to the effective theory, 
Eq. (3.33), and expanding Eq. (3.43) at v = v', one has 

£(1) = 1- (3-44) 

The reduced matrix element £ is the universal function that describes all of the 
matrix elements of operators G between / = | states. It is known as the Isgur- 
Wise function after the discoverers of the relations (3.42) and (3.43). It is quite 
remarkable that the Isgur-Wise function describes both timelike form-factors (as in 
B — > Deu) as well as spacelike form- factors (as in B — > B). The point, of course, is 
that in both cases it describes transitions between infinitely heavy sources at fixed 
"velocity-transfer" (v — v') 2 . 

Expanding Eq. (3.42) for T = 7^ or 7^75, we have 

(D(v')\V^\B(v)) = + t£) (3.45a) 

(DV)e|i M | B(v)) = -^(W)[e*(l + W) -v'^*v} (3.456) 
(D* (v')e\Vfj,\ B(v)) = -£(W)He^ A(J e*^V] (3.45c) 



It remains to express the physical form factors in terms of the Isgur-Wise func- 
tions. We must multiply by \/M D M B to restore to the standard normalization 
of states, and express Eqs. (3.45) in terms of momenta using v = p/M B and 
v' = p'/M D . For example, one has, 

(D(j/)\V v \B(p)) = t{vi/)y/Mrtfo (J?- + j^J (3.46) 



It follows that 



Similarly, /, a± and g can all be written in terms of Moreover, at wo' = 1, 

one has q 2 = (M B v-M D v) 2 = (M B -M D ) 2 = <^ ax so the normalization Eq. (3.44) 
gives 

This remarkable result gives the form factors, in the heavy quark limit, without 
uncertainties from hadronic matrix elements. Short distance corrections will be 
discussed below; see Sections 3.8 and 3.9. 

Exercise 3.7 The same methods can be used to obtain relations among, and normal- 
izations of, the form factors relevant to semileptonic decays of heavy baryons. The case 
of transitions between I = states is simplest. The case of transitions involving higher I 
states can be found elsewhere. [20,17] There are three form factors, Fi, for the matrix ele- 
ment of the vector current between Ab and A c states, and three more, d, for the matrix 
element of the axial current: 

(A c (v', s')\&r„b\A b (v, s)) = « (s ' V)[7„^i + v'„F 2 + v^Fsju^ (v) , 
(A c (v',s')\&y^ 5 b\A b (v,s)) =u is ' ) (v , )[^G 1 +v , tl G 2 + v^G 3 } l5 u is \v) . 

Prove that all six are given in terms of one universal 'Isgur-Wise' function[20]. Show that 
the matrix element of the current is given by 

(A c (v',s')\c v ,Tb v \A b (v, S )) = C(W)« (s V)r« (s) (t>) , (3.49) 
and that ( is fixed at one point: £(1) = 1. Thus (3.49) show that 

Fi = Gi and F 2 = F 3 = G 2 = G 3 = 0, (3.50) 



and Gi(l) = 1 



3.8. Beyond Tree Level 

In the previous sections we have seen how the HQET can be derived from QCD 
and used to obtain useful information about physical processes. The derivation of 
the HQET in Section 3.2 involved only tree level Feynman diagrams. Clearly one 
must extend this beyond tree level if the HQET is to be at all useful. After all, 
we want to use it to describe mesons made of a confined heavy quark and light 
antiquark. It is not difficult to extend the HQET to arbitrary order in perturbation 
theory [6, 21]. We will content ourselves with an understanding of how this works at 
one loop, although going beyond is not much more complicated. 

The alert reader may complain, justifiably, that an all orders proof is not 
enough. There are bona fide non-perturbative effects that one cannot obtain even 
in all orders of perturbation theory. I do not know of a non-perturbative proof 
of the validity of the HQET. One must not forget, though, that many important 
results of quantum field theory are proved perturbatively, e.g., renormalizability of 
the S-matrix in Yang-Mills theories and the operator product expansion. 

The generalization of the effective lagrangian of Eq. (3.1) beyond tree level 
consists of adding to it counterterms, 

^iff = Qviv-DQ V + £i ig ht + £ c .t. • (3.51) 

At tree level the HQET gave us an expression for the Green functions of QCD 
as an expansion of powers in the residual momentum k = p — Mqv. For the two 
point function we had, Eq. (3.11), 

G(p; q) = G v (k; q) + O (A/Mg) (tree level) . (3.52) 

Beyond tree level the corrected version is still close in form to this, 

G(p; q; fi) = C(M q //jl, g s )G v (k; q;fi) + O (A/Mq) (beyond tree level) . (3.53) 

The Green functions G and G v are renormalized, so they depend on a renormaliza- 
tion point fi. The function C is independent of momenta or light quark masses: it 
is independent of the dynamics of the light degrees of freedom. It is there because 
the left hand side has some terms which grow logarithmicaly with the heavy mass, 
ln(Mg//x). The beauty of Eq. (3.53) is that all of the logarithmic dependence on 
the heavy mass factors out. Better yet, since C is dimensionless, it is a function of 



the ratio Mq/h only, and not of Mq and n separately 3 . To find the dependence 
on Mq it suffices to find the dependence on \i. This in turn is dictated by the 
renormalization group equation. 

It is appropriate to think of the HQET as a factorization theorem, stating 
that, in the large Mq limit, the QCD Green functions factorize into a universal 
function of Mq, C(Mq/ '/i, g s ), which depends on the short distance physics only, 
times a function that contains all of the information about long distance physics 
and is independent of Mq, and can be computed as a Green function of the HQET 
lagrangian. 

Exercise 3.8 There is a factorization theorem in the physics of deep inelastic scattering, 
expressing the cross section as a product of parton distributions times parton cross sections. 
Draw an analogy between those quantities in that factorization theorem and the ones in 
the HQET. 

Of course, the theorem holds for any Green functions, and not just for the two 
point function. 

To see how this works at one loop, we start by considering Green functions for 
a heavy quark with n gluons, with n > 1. These are convergent by power counting, 
and since there are no nested divergences at one loop, they are convergent. It 
suffices to consider one-particle irreducible (1PI) functions. In Fig. 1 the left side 
is calculated in the full theory and the right side in the HQET. The double line 
stands for the heavy propagator in the HQET. 




+ 0(1/M) 



Figure 1 Relation between Green functions in full and effective field theories. 

We can prove the validity of the equation represented in Fig. 1, diagram by 
diagram (there are several diagrams that contribute to each side of the equation). 
Consider, for definiteness, the diagrammatic equation in Fig. 2. 



Actually, additional fj, dependence is implicit in the definition of the renormalized 
coupling constant g s . This reflects itself in the explicit form of C. 




+ 0(1/M) 



Figure 2 Relation between Green functions in full and effective field theories in 
simple contribution to four point function, at one loop. 

The equation would trivially hold if we could make the propagator replacement 



even inside the loop integral. Here p = Mqv + k, and / is the loop momentum. In 
other words, in the right hand side of Fig. 2, we take the limit Mq — > oo and then 
integrate, while on the left side we first integrate and then take the limit. Everyone 
knows that, if both integrals converge, then they agree. And that is the case for 
Fig. 2, and, indeed, it is also the case for any 1-loop integral with a heavy quark 
and n > 2 external gluons. We have established Fig. 1 for n > 2. 

We are left with the 2-point (n = 0) and 3-point (n = 1) functions. These are 
different from the n > 2 functions in two ways. First, they receive contributions 
at tree level. And second, they are divergent at 1-loop. Choose some method 
of regularization. Dimensional regularization is particularly useful as it preserves 
gauge invariance (or, more precisely, BRST invariance). The comparison between 
full and effective theories is simplest if the same gauge and regularization choices 
are made. For concreteness, consider Fig. 3. 
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+ 0(1/M) 



Figure 3 Relation between infinite Green functions in full and effective field the- 
ories at one loop: three point function. 



Since both sides are finite, we can argue as before. But we run into trouble 
when we try to remove the regulator. One must renormalize the Green functions by 
adding counter terms, but there is no guarantee that the counterterms satisfy the 
same relation as the regulated Green functions of Fig. 3. To elucidate the relation 
between counterterms, take a derivative from both sides of Fig. 3 with respect to 
either the residual momentum or the gluon external momentum of q^. This 
makes the diagrams finite and the regulator can be removed. Thus, at 1-loop, the 
relations 

" -Gf' 1 ) + O (A/Mq) (3.54) 





d 




dk, 


9 G (2,l) 


d 




dq, 



and t 

" -G^ + 0{A/M Q ) (3.55) 

hold. The counterterms, or at least the difference between them, are k^ and g M 
independent. It is a simple algebraic exercise to show, then, that the difference 
between counterterms is of the form 

aG (2,i)o + 6 £(2,i)o (3.56) 

where the superscript '0' stands for tree level, and a and b are infinite constants, 
i.e., independent of /c M and q^. Thus, one can subtract the 1-loop Green functions 
by standard counterterms, and establish the equality of Fig. 3. 

A similar argument can be constructed for the 2-point function. One must take 
two derivatives with respect to fc M , but that is as it should, since the counterterms 
are linear in momentum. 

We have therefore established that, to 1-loop, the renormalized Green functions 
in the full and effective theories agree. The alert reader must be puzzled as to the 
fate of the function C{Mq/ ji,g 8 ) of Eq. (3.53). What has happened is that the 
constant b in the counterterm in Eq. (3.56) is, in general, Mq dependent. Indeed, 
if we take derivatives with respect to Mq, as in (3.54) or (3.55), the degree of 
divergence is not changed, and one cannot argue that a or b are Mq independent. 
The relation between renormalized Green functions that we have derived contains 
hidden MQ-dependence in the renormalization prescription for the Green functions 
in the HQET. 

Given two different renormalization schemes, the corresponding renormalized 
Green functions G and G' are related by a finite renormalization 



G = z(p,g s )G' 



Choosing G to be the mass-independent subtracted Green function, and G' the 
one in our peculiar subtraction scheme, we have that the relation between full and 
effective theories becomes 

G^\p, q;v)=C (M Q /v, g s ) G^\k, q-») + O (A/Mq) 

as advertised in Section 3.2. Here, C is nothing but this finite renormalization 
z(p,g s ). That we can use the same function C for all Green functions can be 
established by using the same wave functions renormalization prescription for gluons 
in the full and effective theories. Otherwise, an additional factor of z 7 ^ 2 would 
have to be included in the relation between G^ 2,n ^ and G^ 2,n \ This completes the 
argument. 

It is worth mentioning that the discussion above assumes the renormalizability, 
preserving BRST invariance, of the effective theory. Although, to my knowledge, 
this has not been established, there is no obvious reason to doubt that the standard 
techniques apply in this case. 

3.9. External Currents 

We will often be interested in computing Green functions with an insertion of 
a current. Consider, the current 

J r = qTQ (3.57) 

in the full theory, where V is some Dirac matrix, and q a light quark. In the effective 
theory, this is replaced according to 

J r (x) -> e- iM ^ x J r (x) , (3.58) 

where 

J r = qTQ v , (3.59) 

and it is understood that in J r the heavy quark is that of the HQET, satisfying, 
in particular, i/iQ v = Q v . The exponential factor in Eq. (3.58) reminds us to take 
the large momentum out through the current, allowing us to keep the external 
momentum of light quarks and gluons small. The relation between full and effective 
theories takes the form of an approximate equation between Green functions — and 
eventually amplitudes — of insertions of these currents: 

G Jr (p, p'; q- n) = C(M Q /n, g s ) l ' 2 C T {M Q /^ g s )G vJr (k, k'; q; n) + 0(A/M Q ) , 

(3.60) 



where p and p' are the momenta of the heavy quark and the external current, k and 
k' the corresponding residual momenta, p = Mqv + k, p' = Mqv + k' , and q stands 
for the momenta of the light degrees of freedom. The factor C 1 / 2 Cr accounts for 
the logarithmic mass dependence, as explained earlier. We see that an additional 
factor, namely, Cr, is needed in this case to account for the different scaling behavior 
of the currents in the full and effective theories. It is convenient to think of the 
replacement of currents, not as given by Eq. (3.58), but rather by 

J r (x) - e- lM ^ x C T {M Q /^g s )J T {x) . (3.61) 

In fact, Eq. (3.60), and therefore the replacement in Eq. (3.61), are not quite 
correct. To reproduce the matrix elements of the current Jr of Eq. (3.57), it is 
necessary to sum over matrix elements of several different 'currents' in the effective 
theory. The operator Jr of Eq. (3.59) is just one of them. In addition, one may 
have to introduce such operators as qfTQ v . The correct replacement is therefore 

Jr(x) -> e-* M ^J2 C r\M Q /^g s )d^(x) . (3.62) 

i 

Here O^(x) is the collection of the operators of dimension 3 with appropriate 
quantum numbers. The first operator in the sum, call it O (0) , is there even at tree 
level, and corresponds to the operator J r of Eq. (3.59). 

Another case of interest is that of the insertion of a current of two heavy quarks 

Jr = Q'TQ . (3.63) 

The replacement now is 

j r ( x ) ^ e -iM QW +iM g ,»',y-^)^% )W)Jj ^) (l ) (364) 

i V M Q 

Again, O W (x) stands for the complete list of operators of dimension 3 in the effective 
theory with the right quantum numbers. Also, the operator = Q' V ,TQ V appears 
in the sum at tree level. 

This deserves some explanation. The Green functions now include two heavy 
quarks. The functions C connecting these full and effective Green functions will 
now, in general, depend on both Mq and Mq> . Moreover, we can not argue that C 
are independent of the velocities v and v' . In fact, this was true of the simpler case 



considered in Section 3.8; but there, C could only depend on through v 2 = 1. In 
the case at hand there is an additional invariant on which C can depend, namely w) ' . 

The explicit functional dependence on Mq in the functions Cp and Cp can be 
obtained from a study of their dependence on the renormalization point [i. For 
clarity of presentation we neglect operator mixing for now. When necessary, this 
can be incorporated without much difficulty. Taking a derivative d/dfj, on both sides 
of Eqs. (3.53) and (3.61), we find 

/^Cr = (7r-7r)Cr (3.65) 

where 7 r and 7 r are the anomalous dimensions of the currents Jp and Jp in the 
full and effective theories, respectively. Of particular interest are the cases V = 7 M 
and r = 7^75. These correspond, in the full theory, to conserved and partially 
conserved currents, and therefore the corresponding anomalous dimensions vanish, 
giving 

^ = ~ % dr (r = 7 "' r75) • (3 - 66) 

Before we solve this equation, we recall that 

4 = "i +/?w !!' < 3 - 67 > 

Here (3 is the QCD /3-function, with perturbative expansion 

¥--*£ +fc G£0 ,+ - ■ (3 - 68) 

and 

60 = 11-^/, (3.69) 

where rif is the number of quarks in the theory. For our purposes, nj should not 
include the heavy quark. This is explained in the famous paper by Appelquist and 
Carrazone[22]; it simply reflects the fact that the logarithmic scaling of g s is not 
affected by heavy quark loops, since these are suppressed by powers of Mq. Now, 
the solution to (3.66) is standard: 

CrW)=exp(- I" 8 dg'^l) Cp(^9sM) (3.70) 



where g s is the running coupling constant defined by 



» ,d ^f L = mm , 9.W =g.- (3.7i) 

Choosing fi = Mq, and restoring the dependence on Mq, we have then 

C r (M Q /^g s ) = exp f - dg' C r (l,g s (M Q )) . (3.72) 

\ Jg a {M Q ) P{9 ) J 

Therefore, the problem of determining Cf(Mq/ /j, g s ) breaks down into two 
parts. One is the determination of the anomalous dimensions 7 r . The other is the 
calculation of Cr(l, (7 s (Mq)). Both can be done perturbatively, and Cy(Mq/ g s ) 
can thus be computed, provided \x and Mq are large enough so that g s (/x) and 
g s (Mo) are small. One finds, for example, that in leading order Cr is V independent 
and there is no mixing: 

(3 ' 73) 

where a s = g 2 s /Aix, and[23-24] a/ = -c 1 /2b = -6/(33 - 2n f ). 

Exercise 3.9 Obtain an expression for Cr analogous to that of Eq. (3.72) but without 
assuming that the anomalous dimension 7r of Eq. (3.65) vanishes. 

We now turn to the computation of the coefficient Cr for the current of two 
heavy quarks in Eq. (3.64). A new difficulty arises. Because Cr depends on three 
dimensionful quantities, namely the masses Mq and Mq>, and the renormalization 
point /U, its functional dependence is not determined from the renormalization group 
equation (even if we neglect the implicit dependence of g s on ^i). Two different 
approximations have been developed to deal with this problem: 

I) Treat the ratio Mqi/Mq as a dimensionless parameter, and study the de- 
pendence of Cr on Mqi / fj, through the renormalization group[25]. This is just like 
what was done for the heavy-light case, so we can transcribe the result: 

*( i r-s^-)-(-0*) a ( i -s-"'--H- 

(3.74) 

Again 

C r (l,^,tw',&(Mg0) = 1 + 0(a s (M Q/ )) . (3.75) 



But, now, the correction of order 6l s (Mqi) is a function of Mqi/Mq. This method 
has the advantage that the complete functional dependence on Mqi /Mq is retained, 
order by order in 6l s {Mqi). Nevertheless, it fails to re-sum the leading-logs between 
the scales Mqi and Mq, i.e., it does not include the effects of running of the QCD 
coupling constant between Mq and Mqi. Therefore, this method is useful when 
Mqi/Mq ~ 1, or, equivalents, when (a s (M Q i) - a s (M Q ))/a s (M Q ) < 1. 

II) Treat the ratio Mqi/Mq as small. Expand first in a HQET treating Q 
as heavy and Q' as light. The corrections are not just of order A/Mq but also 
Mqi/Mq, but this is assumed to be small (even if much larger than A/Mq). Then 
expand from this HQET, in powers of A/Mq/, by constructing a new HQET where 
both Q and Q' are heavy [16]. The calculation of Cr then proceeds in two steps. 
The first gives a factor just like that of the heavy-light current, in (3.72) 

exp ( - f S(M) dg'^l) C t {1,UMq)) ■ (3.76) 

\ Jg 3 (M Q ) P{9 ) J 

The second factor is as in method I, above, but neglecting Mqi/Mq. Moreover, the 
current Jr is not conserved, so the anomalous dimension to be used is not —% but 
7 r — 7 r . Finally, we must make explicit the fact that in the first and second steps the 
appropriate /3-functions differ in the number of active quarks. We therefore label 
the one in the second step j3' and the corresponding running coupling constant g' s . 
The second factor is 

( /^ s(m) <7 (a') -7 (a")\ - 

exp / dg ~wU - / M'wnA Cr(l,0, W,^ s (Mq,)) • (3.77) 

\Jg s (M Q i) P{9 ) Jg' s {M Q i) P V9 ) J 

Combining factors gives 

Pl [Mqi Mqi . \ ( [ §s{Mq '} ,l r (g') f § '^\ „%(g")\ 

Cr ' ™ >9° ~ exp - / dg' - / dg" j$L± 

\ V Mq J y Jg s (M Q ) P(g') Jg' a {M Q ,) P {9 ) J 

x C r (l,g s (M Q ))C r (l,0, W,g 8 (M Q ,)) . 

(3.78) 

The advantage of method II over method I is that it does include the effects 
of running between Mq and Mqi. The disadvantage is that it neglects powers of 
Mqi/Mq. (Actually, the result can be improved by reincorporating the Mqi/Mq 
dependence, as a power series expansion in this ratio). 



For example, in method II eqn. (3.64) becomes, in leading order, [16] 

c^b - (^\) ai (^rPrY <V[(1 + "h" + (A* - X c (v-v'))^]b v (3.79) 
\a s {m c )J V a s (» / 

for r = 7 M , and 

c 7 ^- (5^)°' f^^) a \v[(l + ^)7^75-(A 6 + A c (W))^ 75 ]6„ 

(3.80) 

for T = 7^75, where 

a s (m b ) , 2a s (m c ) , 

A fo = 3?r , A c (tw ) = — — — r(iw ) , 

g 

a L (v-v') = ^— ^— [v-v'r(v-v') - 1] , (3.81) 

r ( x ) = ; I = ln ( & + \A 2 - l) , 
Vr — 1 v / 

and k is of order a s but a subleading log. 

Exercise 3.10 Why is the term involving k in Eqs. (3.79) and (3.80) a subleading log, 
while those involving and A c are leading logs? 

3.10. Form factors in order a s 

The predicted relations between form factors, and normalizations at g^ ax , are 
only approximate. Indeed, several approximations were made in obtaining those 
results. Corrections that arise from subleading order in the 1/M expansion will be 
considered in Chapter 4. Here we will discuss corrections of order a s . 

As observed in Section 3.9, the vector and axial-vector currents of the full the- 
ory, cTb, match onto a linear combination of 'currents', i.e., dimension 3 operators, 
in the effective theory. At one loop, the correspondence between vector and axial 
currents in the full and effective theories is given by Eqs. (3.79) and (3.80). The 
constant A& and the function A c arise only from 1-loop matching, and are scheme 
independent. The constant k receives contributions both from matching at 1-loop, 
and from 2-loops anomalous dimensions. Leaving out the latter would give a mean- 
ingless, scheme dependent, result. Although k has been computed, it is interesting 
to note that predictions can be made solely form the 1-loop matching computation. 

Indeed, comparing Eqs. (3.45) with Eqs. (3.41), we see that at zeroth order in 
d s (mb) or a s (m c ) we have 

a++a-=0. (3.82) 



Plugging Eq. (3.80) into Eq. (3.42) we see that, to order a s (m c ) and a s (mt,) there 
is a computable correction to this combination of form factors, namely 

a_|_ + a_ ^m c 

a_|_ 777-5 

The constant k, although difficult to compute, does not change the relations be- 
tween form factors since it simply rescales the leading order predictions in Eq. (3.45) 
by the common factor of (1 + k). It does, however, affect the predicted normaliza- 
tion of form factors at q^ax- Since at v' = v the effective vector current is again 
Cv'l^bv, but rescaled by (1 + k + Xb — A c (l)), the correction to Eq. (3.48) is 

We emphasize that retaining the constant k in Eq. (3.84) is inconsistent, because 
not all next to leading logs are included. For a calculation of the sub-leading logs, 
including the constant k, see Ref. [26]. 

4. 1/Mq 

4-1. The Correcting Lagrangian 

One of the main virtues of the HQET is that, in contrast to models of the 
strongly bound hadrons, it lets us study systematically the corrections arising from 
the approximations we have made. To be sure, we've made several approximations 
already, even within the zeroth order expansion in A/Mq. For example, we have 
computed the logarithmic dependence on Mq, i.e., the functions Cp and Cp of 
Eqs. (3.62) and (3.64), using perturbation theory. In this Section we turn to the 
corrections of order A/Mq. 

The HQET lagrangian was derived, in Section 3.2, by putting the heavy quark 
almost on-shell and expanding in powers of the residual momentum, or light 
quark or gluon momentum, q^, over Mq, which we generally wrote as A/Mq. Let 
us again derive the effective lagrangian, keeping track, this time, of the terms of 
order A/Mq. 

We will rederive , including 1/Mq corrections, working directly in configu- 
ration space [27]. The heavy quark equation of motion is 



a s (m b ) 2a s (m c ) , 
H r(w ) 



3tt 



3tt 



(3.83) 



(Hp - Mq)Q = 



(4.1) 



We can put the quark almost on shell by introducing the redefinition 

Q = e~ iM ^ vx Q v (4.2) 

In terms of Q v , the equation of motion is 

[1p + M Q y-l)]Q v = (4.3) 

If we separate the (1 + f) and (1 — f) components of Q v , we see that, as expected, 
the latter is very heavy and decouples in the infinite mass limit. To project out the 
components, 

Qv = Qi +) + Q[- ] (4.4) 

where 



Qi ±} = ( 1 -± l )Q v , (4-5) 



we multiply Eq. (4.3) by {}^jr). Thus we have the equations 

^v■DQ^ = - 0Q<-> (4-6) 

and 



These equations can be solved self-consistently by assuming that Q v is order 
(M Q )° while QV is order Mq 1 . A recursive solution follows. From Eq. (4.7) 

^ , =^(¥)^ +) - i ^- ) < 4 - 8 ' 

Substituting into Eq. (4.6) and dropping terms of order 1/Mq, we have 

w . D Q i+ > = - (l±i) ¥ -L- (izi) ¥ Q i+ > (4.9) 
The right hand side involves 

(4-10) 

where cx^ = §[7^,7^] and G M „ = ^-[.D M ,A,] is the QCD field strength tensor. 
This equation of motion is obtained from the lagrangian 



1 



C^ = Q v iv-DQ v + -^-Q v D 2 -(vD) 2 + ^a^G^ Q v (4.11) 
AMq l z J 

Here I have reverted to the notation Q v for Q v . How to include higher order 
terms in the 1/Mq expansion into should be clear. 



Exercise 4.1 Find the effective lagrangian describing a heavy scalar field to first order 
in 1/M. 



The 1/Mq term in C^ is treated as small. If it is not, it does not make sense to 
talk about a HQET in the first place. It is therefore appropriate to use perturbation 
theory to compute its effects. In this perturbative expansion, the corrections of order 
1/Mq to Green functions, and therefore to physical observables, are computed by 
making a single insertion of the perturbation 

D 2 - {v-Df + |a^Gv] Q v . (4.12) 

The symmetries of the HQET, discussed at length in Sections 3.1 and 3.3, are 
broken by AC. Under the SU(Nf) -flavor symmetry, AC transforms as a combi- 
nation of the Adjoint and Singlet representations, while only the chromo-magnetic 
moment operator 

Qv^G^Q v (4.13) 

breaks the sp'm-SU(2) symmetry: it transforms as a 3 of spin-5't/"(2). 

A single insertion of AC does include all orders in QCD, and it will often prove 
difficult to make precise calculations of 1/Mq effects. Since AC is treated as a simple 
insertion in Green functions, its treatment in the HQET is entirely analogous to 
that of current operators of Section 3.9. There are coefficient functions that connect 
the HQET results with the full theory. It is convenient to include them directly 
into the effective lagrangian as[27,28,29] 

AC = *h-Qv + c 2(vD) 2 + \czg s a^G^] Q v . (4.14) 

ZMq 

Here 

c i = c i {M Q /n,g s ) (4.15) 

can be determined through the methods discussed extensively in Section 3.9. In 
leading- log, one finds [27] 

ci = -1 

8/(33-2?!/) 

" 2 (4.16) 

-9/(33-271/) 



AC 



2 Mr, 



Qv 



c 2 = 3 



C3 = 



a 



■{M Q )J 

M Y 



■{M Q ) J 



4-2. The Corrected Currents 



Just as the lagrangian is corrected in order 1/Mq, any other operator is too. In 
particular, the current operators studied in Section 3.9, are modified in this order. 
At tree level, these corrections are given by the change of variables of last Section: 



J r = qTQ -> qTe 



-iMqvx 



Qv + 



2M Q 



1 - 



tipQv 



(4.17) 



Beyond tree level, this sum of two terms has to be replaced by a more general 
sum over operators of the right dimensions and quantum numbers. The replacement 
is 



7 —IMqvx 



(4.18) 



where Oj are operators of dimension 4 that include, for example, the operators 



qTflpQ v 



qTi(v-D)Q v 



(4.19) 



A complete set of operators, and the corresponding coefficients, Dp , for the cases 
T = 7 M and V = 7^75, can be found in Refs. [28,30] in the leading-log approximation. 
The case of two heavy currents is similar. A straightforward calculation gives 



J r = Q'TQ ^ e -iM Q vx+iM QlV '-x \q' v ,yQ v 



2M, 



Q 



1-4 



2M Q/ 



Q' v >0 



1 - 



(4.20) 



Again, beyond tree level we must replace this expression by a more general sum 
over operators of dimension four, 



j > ^ — iMqvx+iMq/ v'-x 



r °j 



2M Q , 



(4.21) 



It is worth pointing out that, in the computation of the coefficient functions 
D { r j \ D ( r j) and 5 r (j) , there is a contribution from the term of order (1/Mq) . In 
computing the coefficient functions to order 1/Mq one must not forget graphs with 
one insertion of the zeroth order term in the current and one insertion of the first 
order term in the HQET lagrangian. 



4-3. Corrections of order m c /mb 

In the case of semileptonic decays of a beauty hadron to charmed hadron, 
we introduced earlier an approximation method ("Method II" in Section 3.8) in 
which m c /mb was treated as a small parameter. Now, m c /mb ~ 1/3 and you may 
justifiably worry that this is not a good expansion parameter. We will see in this 
Section that the corrections are actually of the order of a s /ir(m c /mb) and therefore 
small. Moreover, they are explicitly calculable. 

The strategy is [28] to look at those corrections of order 1/777,5 which may be 
accompanied by a factor of m c . In the first step of the approximation scheme we 
construct a HQET for the 6-quark, treating the c-quark as light. We must, of course, 
keep terms of order l/m& in this first step. The second step is to go over to a HQET 
in which the m c -quark is also heavy. For now, we care only about terms in this 
HQET that have positive powers of m c . 

In the first step, the hadronic current cT6, with T = 7^ or T = 7^75, is replaced 
according to Eq. (4.18). The question is, which terms in Eq. (4.18) can give factors 
of m c when we replace the c-quark by a HQET quark, c v >. Recall that, once we 
complete the second step, all of the m c dependence is explicit. The answer is that 
any operators in Eq. (4.18) which have a derivative acting on the c-quark will give a 
factor of m c . From Eq. (4.2) we see that a derivative id^ acting on the charm quark 
becomes, in the effective theory, the operation m c v^ + z<9 M . So the prescription is 
simple: take Jr in Eq. (4.18) and replace 

id„ -> m c v'^ (4.22) 

in those terms where id^ is acting on the charm quark. 
For example, if the operator 

— cilpTh 
m h 

is generated at some order in the loop expansion, it gives an operator 

Cy'f 1 V = Cy'L 

after step two is completed. 

It is really interesting to note that the resulting correction does not introduce 
any new unknown form factors. For example, the matrix element of (4.24) between 



(4.23) 



(4.24) 



a B and a D is given by Eq. (3.42) only with an additional factor of —m c /mb in 
front. 

The calculation described here has been performed in the leading-log approxi- 
mation in Ref. [28]. The correction to the vector current is 



AF M = —c v > (ai7 M + a 2 iv + a 3 v' )b v (4.25) 
nib 



where the coefficients = a^/u), written in terms of 



z = c [ (4.26) 
a s (m b ) 



are 



5, , 1 _jl 2W + 12 _jl 34tw'-9 jl 8 , _e_ , 

ai = - W — 1 Z 25 H Z 25 ^25 v . v z 25 i n z 

9 V ; 18 27 54 25 

5,„ rt „ 13 _jl 4W-6 __a. 1W-18 jl 

2 =- 1-2W 2 25 ^ 25 £25 

9 V ' 9 27 27 

15 2 _s_ jl 

a3 — Z 25 — £25 

9 3 



(4.27) 



In particular, this gives a contribution to the form factor, at v = v', of 



vn 

— (a 1 + a 2 + a 3 )U'=i ^- 07 (4.28) 
m & 

This is not negligible! It is reassuring that this type of corrections can be extracted 
explicitly. On the other hand, it should be remembered that both corrections of 
order (m c /m&) 2 and of subleading-log order can still be considerable and should be, 
but have not been, computed. 

4-4- Corrections of order A/ m c andk/rrib- 

Corrections to the form factors for semileptonic decays of S's and A^'s that arise 
from the terms of order l/m c in the effective lagrangian Eq. (4.11) and the currents 
Eqs. (4.18) and (4.21) are, in principle, as large or larger than those considered in the 
previous Section. It is a welcome surprise that the corrections to the combination 
of form factors that contribute to the semileptonic decay vanish at the endpoint 
v-v' = 1. Thus, the predicted normalization of form factors persists, although, as 
we will see, not so the relations between form factors. 

The decay[31] A& — > A c eu is simpler to analyze than the decays[32] B — > Deu 
and B — > D*eu. Moreover, it turns out that for the baryonic decay some relations 



between form factors survive at this order. For these reasons, we will present here 
the baryonic case. We will briefly return to the decay of the meson at the end of 
this section, where we will describe the result. 

There are two types of corrections to consider[31], coming from either the mod- 
ified lagrangian of from the modified current. We start by considering the former. 
The ci and ci terms in the effective lagrangian (4.14) transform trivially under the 
spin symmetry, contributing to the form factors in the same proportion as the lead- 
ing term in Eq. (3.49). This effectively renormalizes the function £ but does not 
affect relations between form factors. 

Moreover, the normalization at the symmetry point wv' = 1 is not affected. 
This is a straightforward application of the Ademollo-Gatto theorem. If is a 
symmetry generating current of a hamiltonian Hq, then corrections to the matrix 
element of the current, at zero momentum, from a symmetry-breaking perturbation 
to the hamiltonian, eHi, are of order e 2 . In the case at hand the Ademollo-Gatto 
theorem implies that corrections to the normalization of £ at the symmetry point 
are of order (l/m c ) 2 . 

The chromomagnetic moment operator in the lagrangian (4.14) does not give a 
contribution at all. The spin symmetries imply 



The function must be an antisymmetric tensor and must therefore be propor- 
tional to v' v v — v^v^. But 



This, we see, is an enormous simplification. There is no analogous simple reason 
for the matrix element of the chromomagnetic moment operator to vanish in the 
case of a meson transition. The chromomagnetic matrix element gives, in that case, 
uncalculable corrections to the relations between form factors. 

We turn next to the contribution from the modification to the current. We 
need the matrix element of the local operators of order l/m c in Eq. (4.21). These 
operators are all constructed of one derivative acting on either heavy quark in the 
quark bilinear. Consider the matrix element 



(A C K, s')\c vl iD^Yb v \K h {v, s)) = u^\v')Yu {s \v)[Av^ + BvJJ, (4.31) 




(4.29) 




(4.30) 



where the form of the right hand side follows again from the spin symmetries. The 
form factors A and B are not independent. Rather, they are given in terms of £. 
To see this, note that, contracting with v'^ and using the equations of motion, 

B = -w'A. (4.32) 

Also, if the mass of the 1 = 1/2 state in the effective theory is A, then 

(A c (v',s')\id^c v ,rb v )\A b (v,s)) = A(« M -t/ Al )C(tw / )«^')(t/)r«W(t;). (4.33) 

Contracting with v^, using the equations of motion and Eq. (4.32) we have 

A(l - (v-v') 2 ) = A(l - vv')C (4.34) 

Therefore, the matrix element of interest is 

(A c (i/, s')\c v nIprb v \A b (r, s)) = A((vv') V " ~ u^(v')-y v Tu^ (v) (4.35) 

1 + v-v 

where V = 7 M or 7^75. Putting it all together, using the leading log expression for 
the coefficients in the current of Eq. (4.21), one finds 

F 1 = Gi r - K < 1 



1 + 

m c V 1 + v-v 

F2 = G2 = - Gl A( 1 ^ (436) 

m c \ 1 + v-v' 

F 3 = G 3 = 



G 1 (D=[^Y (4.37) 



Moreover, 

a s (m c 

as before. Up to an unknown constant, A, there are still five relations among six 
form factors. We can estimate A by writing A = Ma c — m c = (Ma c — M D ) + 
(M D — m c ). If the 'constituent' quark mass in the D meson is ~ 300MeV, then 
A ~ 700MeV. With this, we can estimate the next order corrections to be of the 
order of (A/ 2m c ) -5%. There are, of course, additional computable corrections, 
of order A/2m^ and a s (m c )/n (A/2m c ).[33] 

The result of l/m c corrections to the mesonic transitions is quite different. 
There both the matrix elements of the correction to the current and of the time order 
product with the chromo-magnetic moment operator lead to new form factors. The 
result is that there are incalculable corrections, of order A/2m c , to all the leading 
order relations between form factors. Even if A is smaller in this case, presumably 
A ~ 300MeV, these corrections may be large, say 10%-20%. Remarkably, at the 
symmetry point, v'-v = 1, there are no corrections of order A/2m c to the leading 
order predictions. Thus, one may still extract the mixing angle \V c b\ with high 
precision from measurements at the end of the spectrum of the semileptonic decay 
rates for B — > Deu and B — > D*eu. 



5. Inclusive Semileptonic S-Meson Decays. 

Well before the development of HQETs, it was commonly held that the inclusive 
rate of decay of a heavy meson should be well approximated by the decay rate of 
the heavy quark. Intuitively, the heavy quark just sits at rest (in the rest frame 
of the meson) surrounded by 'brown muck' for which it acts as a static source of 
color. When the heavy quark decays, the complications of strong interactions come 
in, dictating how the rate is divided up between exclusive channels, but the sum 
total must be the quark's rate of decay. 

With the advent of HQET and HQ symmetries one may prove the validity 
of this assertion. Moreover, and more importantly, this is done by setting up an 
expansion for the decay rate in inverse powers of the heavy mass. With a systematic 
expansion one may investigate the accuracy of the result for the not-really-so-heavy 
D and B mesons. 

5.1. Kinematics. 

The decay rate is 

^ = ^£K*w*>'i 2d ** (5.i) 

X 

where j M is the hadronic AB = — 1 current, 

d$ x oc 5^ (P- Px - Pe - Pp ) (5.2) 
is the phase space for final states and 

^ = ^(^)7m( 1 - ls) u {Pe) (5-3) 

is the lepton current. The factor k includes the coupling constants G^|V^| 2 . 
All of the interesting dynamics is in the hadronic tensor 

h^ = J2(^f^\P- Px -q)(m(0)\X}(X\j,(0)\B} (5.4) 
x 

where the sum symbol includes the integral over phase space for the final state X. 

We will investigate the properties of the hadronic tensor h^ u indirectly, by 
studying the related tensor 

Tp V = iJ d%e-^(B(p)\T(jt(x)UO))\B) (5.5) 



Both tensors can be written in terms of Lorentz invariant form factors, thus: 

h^u = guvhi + p^p u h 2 H (5.6a) 

T^u = g^Ti + p^p v T 2 H (5.66) 

Exercise 5.1 I have intentionally left the list of form factors incomplete. The ellipsis in 
Eqs. (5.6) stand for a finite number of terms. Ennumerate them. There is one term that 
can be eliminated by the requirement of time reversal invariance of the strong interactions. 
Which one is it? 

The form factors Tj are more amenable to theoretical investigations because 
is a Green function. The hi form factors can be obtained from these by cutting: 

hi = -ImTj for p ■ q < m| (5.7) 

TT 

Exercise 5.2 To see this, 

(i) Separate the integral in Eq. (5.5) into a; > and x° < parts. Insert a complete set 
of states between the two currents, and make the dependence on x explicit. 

(ii) Carry out the integrations, using, eg, 



/ 



d A x 6(x°y (p - px - q) x = (2tt) 3 5(P -px-q)- 



P°-p° x -q° + ie 
and use the standard identity 



z ± it z 

where PP stands for the principal part, to separate "imaginary" from "real" parts 
(these are in quotation marks because by them I just mean the parts that come from 
the PP terms and from the 5 term, respectively). 

(iii) In the "imaginary" part you will recognize one term is just h^w Write out explicitly 
the other term in the "imaginary" part and show that it vanishes if P ■ q < vtlb- (Hint: 
What is the lowest mass of the intermediate state X?). 

(iv) Complete the proof by relating form factors and taking, where appropriate, the real 
or imaginary parts. 



5.2. The Analytic Structure of The Hadronic Green Function 

The hadronic form factors Tj and hi are Lorentz invariant functions of the 
momenta P and q. There are only two invariant variables, namely Q = a/q 2 /m 2 B 
and z = P • q/m 2 B . Here and throughout a hat on a dimensionful quantity means 
taking out the dimensions by apropriate powers of tub- We will study the behaviour 
of T^ for fixed Q 2 as a function of z. This variable is of interest because 1 — z 
measures how far one is from the resonant region. To see this, consider any final 
state X contributing to the inclusive decay, so q = P—px- Then z = 1 — P-px/rn^, 
or, in the rest frame of the decaying meson, 1 — z = Ex/rriB- Only states with 
small energy contribute to the region of z ~ 1 , and only single resonances contribute 
when z is large enough. 



Exercise 5.3 Show that the physical region extends from z = (neglect the lepton mass) 
to 

z = |(1 + Q 2 - ml) = zq 



Since T Mi , is a Green function we can study its behaviour as a function of 
complex momenta. In particular we can hold Q fixed and study T{(z) as functions 
of complex z. What do we know about the complex z plane? You established above 
that the segment of the real line (0, zq) corresponds to the physical region for the 
inclusive semileptonic B decay. And we know that the right end of the segment 
corresponds to the resonance region. See Fig. 4. 

Im z 



Resonance Region 




Rez 



Physical Region 



Figure 4 The location of the physical and resonant regions in the complex z-plane. 
The variable z is the energy of the final hadronic system in the semileptonic B 
decay, in the B restframe. 

Associated with the physical process there is a cut in the z plane. The cut 
extends from zq = | (1 + Q 2 — rh 2 ) to — oo on the real axis. There is, in addition, a 
cut that extends from z\ = -|(3 — Q 2 ) out to +oo on thereal axis. See Fig. 5. 

Exercise 5.4 What are the physical processes that correspond to the cuts in the z plane? 

If the presence and location of these cuts does not seem obvious, the reader may 
use Landau conditions 4 to determine the analytic properties of a Green function in 



4 See, for example, Ref. [18] 



perturbation theory. She will then find the same cuts, save for the fact that the 
masses will turn out to be those of quarks rather than of the physical bound states. 



Figure 5 The location of the cuts in the complex z-plane. 

As will become evident later, a perturbative computation will be trustworthy 
provided z is nowhere near the resonance region. We should stay away from the 
resonance region by at least a typical hadronic scale, say the rho mass, m p , or the 
chiral symmetry breaking scale, A x . This is why the problem at hand is non-trivial: 
to compute the form factors hi from the imaginary part of a Green function one 
needs to deal with the resonant region where perturbation theory is not valid. 

The situation is much better if we are willing to settle for slightly less, namely 
for averages over the variable z of the decay rate (at fixed Q). Integrating both 
sides of Eq. (5.7) one has 



Im z 



Rez 




(5.8) 



where the contour C runs above the cut from to zq, around the branch point and 
back to under the cut; see Fig. 6. 
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Figure 6 Contour of integration, C, in the dispersion relation (5.8). 

We can trade the contour C for a new one C that stays away from the resonance 
region. C is restricted to start just above the origin of the complex z plane and to 
finish just under it. This is justified because there are no singularities other than the 
cuts. See Fig 7. Provided zq and z\ are well separated there is no problem getting 
the contour C to lie well away from the cut, exept in the vicinity of z = 0. But in 
particular one can have the cut away from the resonance region in its entirety, and 
one may apply perturbation theory to the computation of the Green function. 

Exercise 5.5 Under what circumstances can the cuts meet? In other words, determine 
the kinematic regime for which z\ < zq. What does this imply for inclusive b — > cev and 

b — > uevl 

The only thing missing is a means of computing the Green function perturba- 
tively. The problem is that we need to compute a matrix element between meson 
states, and this usually involves complicated dynamics. As we will see in the next 
section it is possible to use the HQET and the HQ symmetries to compute the 
matrix element. Before turning to that problem, we conclude with a simple obser- 
vation that extends the above result. The average in Eq. (5.8) uses one of many 
possible measures, namely simply dz. One could just as well consider 



I dz f(z)hi(z) 
Jo 



(5.9) 



Clearly we want the function f(z) to be regular on the region of integration. More 
generally we would like its singularities in the complex z plane to be as far away 
from the contour C as possible. This suggest having them at infinity, ie, using an 
entire function. Particularly interesting is the set of functions f n {z) = (zq — z) n , 
which measure the relative importance of the non-resonant and resonant regions. 
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Figure 7 Deformed contour of integration, C, in the dispersion relation (5.8). 
Note that it stays away from the physical region everywhere except near z—0, 
provided Zq and Z\ are well separated. 

5.3. An HQET based OPE 

The Operator Product Expansion (OPE) has been used with great success 
to separate perturbative from non-perturbative effects in deep inelastic nucleon 
scattering. The coefficient functions may be computed perturbatively, while the 
non-perturbative information is encoded in the matrix elements of the operators. 

The common OPE cannot be applied to T^ u . The problem is that it is an 
expansion in inverse powers of Q. The matrix elements of operators of heavy mesons 
can in fact grow with powers of the heavy mass tub- Since in the physical region 
Q < rriB, the ratio uib/Q is not a good expansion parameter. 

But there is a way one can produce a useful OPE for T^ u . The idea is simple: 
expand in inverse powers of the largest scale around, namely mg. Clearly we should 
be using the techniques of the HQET to produce such an expansion. Let P = msv. 



Consider any matrix element (physical or unphysical) of the time ordered product 
T(jl(x)j fl (0)). The idea is to take the momentum of the b quark to be p = m b v + k 
and to expand in inverse powers of m&. 

Consider, in particular, the matrix element between b quarks of momentum p. 
The tree level Feynman graph is in Fig. 8. We insist that the momentum in the 
intermediate light quark, q = u or c, be m^v + ■ ■ ■ and expand in inverse powers 
of the mass. Arbitrary quantities, like Q 2 , may scale with m b (so we hold, eg, Q 2 
fixed as m& — > oo). 

u,c 





mv+k mv+k 

Figure 8 Tree level Green function for the product of currents. 

To see how the expansion works, let us compute the Feynman diagram of Fig. 8. 
To be concise we consider the current jr = qTb, where V stands for any Dirac gamma 
matrix. Then we have the tree level Green function and its expansion: 



1 r- 



m h j + $ - 4 - m q m b j - 4 - m q + # /m b 



m b i)-i-rh q 



$ 1 ^ 



(5.10) 



This matrix element can then be used to determine the coefficients of operators in 
an OPE-like relation between operators: 

d A x e- iq - x T(bYq(x) qTb(O)) 



= — i v r 1 - rb v - —l v r 1 — up - 1 - r&„ + • • 

mb I f - q- m q m b f - q - m q p-q-m q 

(5.11) 

where the ellipsis stand for an expansion in powers of l/m b which has operators 
of ever increasing dimensions. Going beyond tree level will modify the specific 
coefficient functions but the form of the expansion will remain practically the same 
(some new operators may be introduced beyond tree level). 




mv+k mv+k 

Figure 9 Representation of the combined OPE-HQET at tree level. 

Note that the right hand side of Eq. (5.11) is written in terms of the quark 
operators of the HQET, while the left hand side involves the original quark fields. 
This is what makes this OPE different from the common one (used, say, in deep 
inelastic scattering). And there is no way out of it: since we are expanding the full 
theory time ordered product in powers of l/nib, we must use HQET fields for the 
result. A useful pictorial representation of this is presented in Fig. 9, which shows 
the HQET-OPE at tree level. The diagrams on the right hand side of the equation 
have double lines, reminding us they represent quarks in the HQET. 

The coefficients in the OPE of Eq. (5.11) have poles in z — zq of increasingly 
higher order. Symbolically we may write for the expansion of the operator product 

m b ^ {2m b ) n 
with coefficient functions that behave as 

(5.13) 



(z - z 



\n+l 



Beyond tree level there will be cuts in additions to these poles. Although the nature 
of the singularities change, the discussion that follows remains essentially the same. 

The first two terms in the OPE expansion of T Mi , can be computed! The first 
one involves the matrix element between S-mesons of a dimension-3 operator. But 
this is fixed by spin/flavor symmetries: 

(B\b v T^+4 + m q )Tb v \B) = -Z(l)tTB(v)ry + 4 + m q )rB(v) (5.14) 

with £(y-v') the Isgur-Wise function, £(1) = 1. The second one involves a dimension- 
4 operator. The matrix element of this vanishes: 

(B\b v TD^b v \B) = (5.15) 



There are therefore no corrections of order l/nib- More on this below. 



Exercise 5.6 Prove Eq. (5.15). Hint: Consult section 4.4. 



We see that one can compute the rate for inclusive semileptonic B decay. With 
the tree level coefficient functions just computed, the matrix elements in Eqs. (5.14) 
and (5.15), and neglecting higher order terms in l/m&, one obtains for the decay 
rate the same result as if computing the decay rate of a free b quark! 

Exercise 5.7 Prove this assertion. In fact, we can do even better because we need not 
integrate over Q 2 . Prove that the moments 



are the same as for the free quark decay. This is an important exercise. You will have 
to put together all the pieces of the argument discussed above, and then more. 

The deviations from the free decay rate are parameterized by the matrix ele- 
ments of the operators of dimension 5 and higher. In this sense we have managed 
to separate long from short distance effects. 

Note that for z on the real axis, the expansion breaks down as z — > z . More 
precisely, in this limit the coefficients c n of Eq. (5.12) grow arbitrarily large. This 
is in accord with the statement we had made previously that the perturbative 
calculation of the Green function breaks down in the resonant region, z ~ zq. 

A physically interesting choice of the measure / in Eq. (5.9) is /i = z — zq, 
because it measures the deviation from the free quark result. The quantity m\f\(v ■ 
q/rrib) is the invariant mass-squared of the hadronic final state minus m 2 q . If the 
final state quark produced by the b decay is on its mass shell, then fi = 0. 

Exercise 5.8 Prove that with fi = z — zq the contribution from the parton model to 
the average in Eq. (5.9) vanishes. Discuss corrections to this result arising from the 
nonperturbative terms and from the perturbative corrections to the coefficients in the 
expansion of Eq. (5.12). 



Exercise 5.9 \V u b\ is determined from a measurement of inclusive charmless semileptonic 
B decay. Since V u b\ ^ \ Vcb\, the rate for B — > Xei> is dominated by the decay into charm. 
Short of reconstructing the hadronic state X, one is forced to consider the spectrum of 
electron energies, E e . Find the maximum electron energies, E^ ,xaax an d i?e' max , in the 
decay into charm and charmless final states, respectively. What does the region E^' max < 
E e < Eg ,max , correspond to in terms of our variables z and Q? 




6. Chiral Lagrangian with Heavy Mesons 



6.1. Generalities 

Chiral symmetry and soft pion theorems have been used in particle physics for 
several decades now with great success. The most efficient way of extracting infor- 
mation from chiral symmetry is by writing a phenomenological lagrangian for pions 
that incorporates both the explicitly realized vector symmetry and the non-linearly 
realized spontaneously broken axial symmetry. [34] Theorems that simultaneously 
use heavy quark symmetries and chiral symmetries are most expediently written by 
means of a phenomenological lagrangian for pions and heavy mesons that incorpo- 
rates these symmetries. [35,36] 

In the limit m b — > oo, the B and the B mesons are degenerate, and to imple- 
ment the heavy quark symmetries it is convenient to assemble them into a "super- 
field" H a (v): 



Ha(v) = -^-[B a ^-5a7 5 j • (6-1) 

Here is the fixed four-velocity of the heavy meson, and a is a flavor SU(3) 
index corresponding to the light antiquark. Because we have absorbed mass factors 



\j2m,B into the fields, they have dimension 3/2; to recover the correct relativistic 
normalization, we will multiply amplitudes by for each external B or B 

meson. 

The chiral lagrangian contains both heavy meson superfields and pseudogold- 
stone bosons, coupled together in an SU(3) l x SU(3)r invariant way. The matrix of 
pseudogoldstone bosons appears in the usual exponentiated form £ = exp(iA^//), 
where 

-73*° + > K ° I, (6.2) 



M 



n 



■o 



and / is the meson decay constant. The bosons couple to the heavy fields through 
the covariant derivative and axial vector field, 



ab 



K b = \ - ^) ab = -jd,M ab + 0(M 3 ) • 



Lower case roman indices correspond to flavor 577(3). Under chiral SU(3)l x 
SU(3)r, the pseudogoldstone bosons and heavy meson fields transform as 



H -> i7[7 f 
(TJ^/7) -> (D^H)U^ 

where the matrix £7 is a nonlinear function of the pseudogoldstone boson matrix 
«M, implicitly defined by the transformation law for £. Then E = ^transforms as 
follows: 

E -> LEi? f . 

The effective lagrangian is an expansion in derivatives and in inverse powers of 
the heavy quark mass. The kinetic energy terms take the form 

£ kin = \f 2 d^ ah d^\ a - Tr \H a {v)iv ■ D ba H b (v)} . 

Here the trace is in the space of 4 x 4 Dirac matrices that define the "superfields" 
H a (v) in Eq. (6.1). The leading interaction term is of dimension four, 

£ int = gTr \H a (v)H b (v) 4 6a75 ] , (6.3) 

where g is an unknown parameter, of order one in the constituent quark model. 
The analogous term in the charm system is responsible for the decay D* — > Dtt. 
Expanding the term in the lagrangian in (6.3) to linear order in the Goldstone 
Boson fields, M, we find the explicit forms for the D*DM and D*D*M couplings 

D*»d,Mrf + h.c. + (^j e^ XK D*»d»MD* x v K . (6.4) 

Using this one can compute the partial width 

r(D*+ -. D+n°) = 

The ACCMOR collaboration has reported an upper limit of 131 KeV on the 
D* width. [37] The branching fractions for D* + — > 7J°7r+ and D* + — > 7J+7T are 



(68.1 ± 1.0 ± 1.3)% and (30.8 ± 0.4 ± 0.8)%, respectively, as measured by the CLEO 
collaboration. [38] Using / = 130 MeV, one obtains the limit g 2 < 0.5. Even if the 
D* decay width is too small to measure, radiative D* decays provide an indirect 
means for determining the coupling g, and provide a lower bound (7 2 >0.1.[39] 

Since charmed and beauty baryons are long lived, one can write down phe- 
nomenological lagrangians for their interactions with pions. These are as well jus- 
tified and should be as good an approximation as the lagrangian for heavy mesons 
discussed above. The treatment is rather similar, and we refer the interested reader 
to the literature. [40] 
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Figure 10 Feynman diagrams for B — > Deu 
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6.2. B -> Deu And B -> D*nEu 

As a first example of an application consider a soft pion theorem that relates 
the amplitudes for B — > D*eu and B — > D*neu, for small pion momentum. [36] The 
heavy quark current is represented in the phenomenological lagrangian approach by 

4 b = - To)h^ - -£(W)TrFiV)7M(i - To )Hi b \v) + ■■■ (6.5) 

where the ellipsis denote terms with derivatives, factors of light quark masses m q , 
or factors of 1/Mq, and £(W) is the Isgur-Wise function. 5 The leading term in 

5 The symbol £ is traditionally used both for the Isgur-Wise function and for the expo- 
nential of the meson fields. To distinguish between them, whenever context may not be 
sufficient, we denote the former as a function of velocities, £(m/). 



Eq. (6.5) is independent of the pion field. Therefore, it is pole diagrams that 
dominate the amplitude for semileptonic B — > Dn and B — > D*n transitions; see 
Fig. 10. These pole diagrams are calculable in this approach, and are determined 
by the Isgur-Wise function and the coupling g. 
A straightforward calculation gives [36] 

{D{v')-K a {q)\jf\B{v)) = iu{BY\T a u{D)^M^M^^{v-v') 

x I — [ie^\ K q"v ,x v K + q ■ (v + v')v^ - (1 + 

I v ~q 

-^[ie^\ K q J/ v' x v K + q ■ (v + - (1 + 

where u(M) stands for the isospin wavefunction of meson M. A similar but lengthier 
expression is found for B — > D*rreu. Even if the coupling g is close to its upper 
limit, this expression makes a small contribution to the inclusive semileptonic rate. 

6.3. Violations To Chiral Symmetry 

Phenomenological lagrangians are particularly well suited to explore devia- 
tions from symmetry predictions. We begin by introducing symmetry break- 
ing terms into the phenomenological lagrangian. The light quark mass matrix 
m q = diag(m u , rrid, m 8 ) parametrizes the violations to flavor SU(3)v- To linear 
order in m q and lowest order in the derivative expansion, the correction to the 
phenomenological lagrangian is 

AC = A [rriqY, + m g S t ] a 

+ X 1 Tr H iQ)a H^ [tm^ + ^m^fa 

+ X' 1 TrH (Q)a HP [m q E + m q ^] b b (6.6) 

The coefficients Ao, Ai and A^ are fixed by non-perturbative strong interaction 
effects, but may be determined phenomenologically. We postpone consideration of 
mass relations obtained from this lagrangian until we have introduced heavy quark 
spin symmetry breaking terms into the lagrangian too. 

Exercise 6.1 The operator \^m q £,—^m q ^^ is parity odd. One could add a term 
TrF (Q)a ^ Q) 7s [£m q £ - ^m q ^] b a . Why did I leave it out? 



The decay constants for the D and D s mesons, defined by 

(0\d^ l5 c\D+(p))=if D p^ (6.7) 

and 

(0\s~/^ 5 c\D s (p)) = if Ds Pn , (6.8) 

determine the rate for the purely leptonic decays D + — > [i^v^ and D s — > /U + z/ M . 
These are likely to be measured in the future. [41] In the chiral limit, where the up, 
down and strange quark masses go to zero, flavor SU(3)v is an exact symmetry 
and so fo s /fD+ = 1- However m s ^ 0, so this ratio will deviate from unity. 

In the chiral lagrangian approach the current operators in Eqs. (6.7) and (6.8) 
are represented by 

Am) = y ^[7 A 75^)4] (6-9) 

This is the lowest order only in an expansion in derivatives. Chiral symmetry vio- 
lation in the relation between fo s and fjy+ arises directly from symmetry breaking 
terms that can be added to the current in Eq. (6.9) and indirectly in loop diagrams 
through the symmetry breaking in the masses of the virtual particles. The latter 
involves, at one loop, the Feynman diagrams in Fig. 11, where a dashed line stands 
for a light pseudoscalar propagator. Neglecting the up and down quark masses in 
comparison with the strange quark mass, the loop graphs give[42,43] 

fD s /f D+ = 1 - | (1 + 3s 2 ) In (Ml/^) +■■■ (6.10) 

Here fi is the renormalization point. The contribution from r\ loops has been written 
in terms of Mk using the Gell-Mann-Okubo formula = 4M| r /3, and the contri- 
bution from pion loops, proportional to lnM^, has been neglected. The ellipsis 
denote terms proportional to the strange quark mass (recall M\ ~ m s ) without a 
log. Theoretically, as m s — > the log term, which has non-analytic dependence on 
m 8 , dominates. The neglected terms, as well as the contribution from the symmetry 
breaking terms in the current are analytic in m s . 

Exercise 6.2 Display the symmetry breaking additions to the current in Eq. (6.10) to low- 
est order in m q . This will introduce new unknown parameters. Compute the modification 
to /d s / fo+ ~ 1 from these terms. 



The dependence of the new parameters in the current on the subtraction point 
\i cancels that of the logarithm. Without additional information on these new 
parameters we have no real predictive power. But one can venture a guess by 
arguing, as above, that the log term dominates. What value should we assign the 
renormalization point \xl If \x is of order the chiral symmetry breaking scale then the 
new parameters have no implicit large logarithms. Numerically, using \i = 1 GeV, 
the result is that 

f Ds /f D+ = 1 + 0.064 (l + 3# 2 ), (6.11) 
or f D Jf D + = 1-16 for g 2 = 0.5. 
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Figure 1 1 Feynman Diagrams in the calculation of Jd s / fp . 

The same formula also holds for fsJfB- In fact, to leading order in 1/Mq the 
ratio is independent of the the flavor of the heavy quark. Consequently, 

IbJIb = 1 (6J2) 

JdJJd 

to leading order in 1/Mq and all orders in the light quark masses. Now, Eq. (6.12) 
also holds as a result of chiral symmetry, for any m c and mj. That is JbJ Jb 
and /d s //d are separately unity in the limit in which the light quark masses are 
equal. This means that deviations from unity in Eq. (6.12) must be small, 0(m s ) x 
0(l/m c — l/m&).[44] This ratio of ratios is observed to be very close to unity in a 
variety of calculations. [45] This may be very useful, since it suggests obtaining the 
ratio /s s //b of interest in the analysis of B — B mixing (see below) from the ratio 
/d s //d, measurable from leptonic D and D s decays. 

Exercise 6.3 The hadronic matrix elements needed for the analysis of B — B mixing are 
(B(v)\b^(l - j 5 )d b^(l - 75)d|5(«)> = p 2 B B B , 

(5 a (v)\h"(l ~ I5)s 67^(1 " 75)s\B s (v)) = p 2 Ba B Bs , 

where the right hand side of these equations define the parameters B Ba and B B . In the 
SU(3)v symmetry limit B Ba /B B = 1. For non-zero strange quark mass, the ratio is no 
longer unity. Show that the chiral log is 

^-l-§(l-*-)^l»(*i/M") + - 

Again, M 2 = 4M^/3 has been used. Using u, = 1 GeV, / = fx, and g 2 = 0.5, the 
correction is B Bs /B B ~ 0.95. Hint: The hard part is to represent the four-fermion 
operators in the effective theory. If stuck, consult Ref. [42] 



Violations to chiral symmetry in B — > D semileptonic decays have also been 
studied. One obtains that a different Isgur-Wise function must be used for each 
flavor of light spectator quark [43] 

igl = 1 + l9 2 ^')j^p I" (Mllf) + *(/•. W)Ml + ■ • . (6.13) 
where A'(/i, m/) is the analytic counter-term, and 



. 2 + x / a; + 1 + \/a; 2 — 1 \ a; , / x — y/x' 

Q(x) = -1 + In + In 



2 V^ 2 — 1 \ x + 1 — Va; 2 — 1 / 4 Va^ 2 — 1 \ a: + v^c' 
or, expanding about a: = 1, 

0( X ) = -I (a; _ i) + —( X - l) 2 - —( X - l) 3 + . . . 
v ' 3 V ; 15 V ; 35 V ' 

Using g 2 = 0.5 and \i = 1 GeV, and neglecting the counterterm one obtains 

Uv-v') 



fu,d(w>') 

or a 5% correction at w/ = 2. 



1 -0.21O(W) + 



Exercise 6.4 Display the Feynman graphs that contribute to the ratio in Eq. (6.13). 

6.4- Violations To Heavy Quark Symmetry 

In a similar spirit one can consider the corrections in chiral perturbation theory 
to predictions that follow from heavy quark spin and flavor symmetries. These are 
effects that enter at order 1/Mq, so the first step towards this end is to supplement 
the phenomenological lagrangian with such terms. In particular, there are SU(3)v 
preserving terms of order 1/Mq that violate spin symmetry in the lagrangian, such 
as [42, 46] 

A£ int = -^-TrH iQ)a a^HPa^ + j^-Tr \H a {v)A hal ,H h {v)] . (6.14) 

In addition there are contributions to the lagrangian in order 1/Mq that violate 
flavor but not spin symmetries. These can be characterized as introducing Mq 
dependence in the couplings g, X± and A^ of Eqs. (6.3) and (6.6). At the same order 
as these corrections, there is a term that violates both spin and SU (3)y symmetries 

A 3 



'^V^v] fem^ + eW*] 6 . (6-15) 



Exercise 6.5 Make a complete analysis of terms that violate SU(3)v and heavy-quark 
symmetries simultaneously, to leading order in both m q and 1/Mq. Note that our list is 
shorter since we are going to concentrate shortly on the spectrum only. 



Spin symmetry violation is responsible for "hyperfine" splittings in spin mul- 
tiplets. To leading order these mass splittings are computed in terms of the spin 
symmetry violating coupling of Eq. (6.14) 

Q \ 

A B = M B * -M b = (6.16) 

m b 

That the mass splittings scale like 1/Mq seems to be well verified in nature: 

M D * - M D M B 
M B * - M B ~ M D 



Measured Mass Splittings 



X - 


Y 


M x -My(MeV) 


D s - 


-D+ 


99.5 ± 0.6[47] 


D+ 


-D° 


4.80 ± 0.10 ±0.06[48] 


D*+ 


-D*° 


3.32 ± 0.08 ±0.05[48] 


D*° 


-D° 


142.12 ± 0.05 ±0.05[48] 


D*+ - 


-D+ 


140.64 ± 0.08 ±0.06[48] 


D* 8 - 


-D s 


141.5 ± 1.9[47] 


B s 


-B 


82.5±2.5[48] or 121 ±9[49] 


B°- 


-B+ 


0.01 ± 0.08[47] 


B* 


-B 


46.2 ± 0.3 ± 0.8[50] or45.4 ± 1.0 [49] 


Bt- 


-B s 


47.0 ± 2.6[49] 


(D*° 


-D°) 




-(D*+ 


-D+) 


1.48 ± 0.09 ±0.05[48] 



Armed with the machinery of chiral lagrangians that include both spin and chi- 
ral symmetry violating terms, one can compare hyperfine splitting for different fla- 
vored mesons. There is a wealth of experimental information to draw from; see Ta- 
ble 3. Breaking of flavor 577(3) y and heavy quark flavor symmetries by electromag- 
netic effects is not negligible. It is readily incorporated into the lagrangian in terms 
of the charge matrices Qq = diag(2/3, —1/3) and Q q = diag(2/3, —1/3, —1/3), [51] 
which must come in bilinearly. For example, terms involving Q 2 q correspond to re- 
placing m q — > Q q in Eqs. (6.6) and (6.15). The electromagnetic effects of the light 
quarks can be neglected if one considers only mesons with d and s light quarks. The 



d ) 



electromagnetic shifts in the hyperfine splittings A Xq and A Xq (X = D , B , q = d, s) 
differ on account of different b and c charges, but they cancel in the difference of 
splittings 

A Xs - A Xd = (M x . -M Xs )- (M x * - M Xd ) 

The only term in the phenomenological lagrangian that enters this difference is 
Eq. (6.15). This immediately leads to 

_9_ 

(Mb* -M Ba )- (M B *-M Bd ) = ^ (^fyj [(Md; - M D .) - (M D * - M D 

(6.17) 

We have included here the short distance QCD effect. [27] 

One expects Eq. (6.17) holds to much better accuracy than the separate rela- 
tions for each hyperfine splitting in Eq. (6.16). Recall that SU(3)v breaking by 
light quark masses and electromagnetic interactions have been accounted for in 
leading order. Moreover, the result is trivially generalized by replacing the quark 
mass matrix in Eqs. (6.6) and (6.15), by an arbitrary function of the light quark 
mass matrix. It is seen from Table 3 that this relation works well. The left side is 
1.2 ± 2.7 MeV while the right side is 3.0 ± 6.3 MeV. 

Since both sides of Eq. (6.17) are consistent with zero and both are propor- 
tional to the interaction term in Eq. (6.15), it must be that the coupling A3 is very 
small. [51] From the difference of hyperfine splittings in the charm sector 

sX3 (m s - m d ) = 0.9 ± 1.9 MeV 



while 

12A 

M D - M Dd = 4Ai(m s - m d ) -(m s - m d ) = 99.5 ± 0.6 MeV 

m c 

leading to | A3/A1 j less than ~ 20 MeV. This is smaller than expected by about an 
order of magnitude. With such a small coefficient it is clear that the next-to-leading 
terms and the loop corrections may play an important role. In particular they may 
invalidate the simple 1/Mq scaling of Eq. (6.17). [52] There is no obvious breakdown 
of chiral perturbation theory, even though the leading coupling (A3) is anomalously 
small. [53] 

At one loop, the expressions for the mass shifts involve large 0(m s lnm s ) and 

3 /2 

0(m s ) (non-analytic) terms. [43,53] The coupling A3 is not anomalously small at 
one loop. Instead, the smallness of the difference of hyperfine splittings in Eq. (6.17) 



is the result of a precise cancellation between one loop and tree level graphs. Ex- 
plicitly, [53] 



With g 2 = 0.5 and \x = 1 GeV, the chiral log is 30 MeV, so the A3 counterterm 
must cancel this to a precision of better than 10%. 

The 1/Mq corrections to the masses Mx and Mx* drop out of the combination 
M x + 3M X *. The combination (M Xs + 3M X *) - (M Xd + 3M X *) is a measure of 
SU(3)v breaking by a non- vanishing m s (or m s — if the d quark mass is not 
neglected). It can be computed in the phenomeno logical lagrangian. To one loop [53] 



\ (Mx s + 3M X ;) - \ (M Xd + 3M X .) = 4A!m s - g 2 (l + JLi) 

(6.18) 

The pseudoscalar splittings (Mo a — Mo d ) and (Mb 3 — MB d ) have been measured; 
see Table 3. Also, \ (M Xs + ZM X . )-\(M Xd + 3M X . ) = f [{M x . -M Xa )- (M x * - 
Mx d )] + {Mx a — Mx d ), and the term in square brackets is less than a few MeV, as 
we saw above. The combination (Mx s + 3Mx*) — (Mx d + 3Mj«) in Eq. (6.18) is 
first order in m s but has no corrections at order 1/Mq. Thus, one expects a similar 
numerical result for B and D systems. Experimentally, (Mb s —MB d )/ (Md s —Mo d ) 
is consistent with unity; see Table 3. The formula in Eq. (6.18) has a significant 
contribution from the term which is independent of the splitting parameter 
Ai. The M\ term gives a negative contribution to the splitting of ~ —250 MeV 
for g 2 = 0.5. The chiral logarithmic correction effectively corrects the tree level 
value of the parameter Ai; for ji = 1 GeV and g 2 = 0.5, the term 4Aim s gets a 
correction 0.9 times its tree level value. Thus, the one-loop value of 4Aim s can 
be significantly greater than the value determined at tree-level of approximately 
100 MeV. 

Chiral perturbation theory can be used to estimate the leading corrections to 
the form factors for semileptonic B — > D or D* decays which are generated at low 
momentum, below the chiral symmetry breaking scale. Of particular interest are 
corrections to the predicted normalization of form factors at zero recoil, v-v' = 1. 
According to Luke's theorem (see Section 4.4), long distance corrections enter first 
at order 1/Mq. Deviations from the predicted normalization of form factors that 



arise from terms of order 1/Mq in either the lagrangian or the current are dictated 
by non-perturbative physics. But there are computable corrections that arise from 
the terms of order 1/Mq in the lagrangian. These must enter at one-loop, since 
Luke's theorem prevents them at tree level, and result from the spin and flavor 
symmetry breaking in the hyperfme splittings and Ag. Retaining only the 
dependence on the larger Ad, the correction to the matrix elements at zero recoil 
are [54] 

(D(v)\J*\B(v)) 

(D*(v,e)\J*\B(v)) 

(6.196) 

where C and C stand for tree level counter-terms and 

f°° z 4 ( 1 1 \ 

F( ~ X) ^J ^(^2 + 1)3/2 \[( z 2 + 1) l/2 + x] 2 -^TlJ (6,20) 

As before, no large logarithms will appear in the functions C and C if one takes 
\i ps Ail f ~ 1 GeV. With this choice, formally, their contributions are dwarfed by the 
term that is enhanced by a logarithm of the pion mass. Numerically, with g 2 = 0.5 
the logarithmically enhanced term is —2.1% and —0.7% for D and D* , respectively. 

The function F accounts for effects of order (l/m c ) 2+n , n = 1,2, ... It is 
enhanced by powers of 1/M n over terms that have been neglected. Consequently it is 
expected to be a good estimate of higher order l/m c corrections. With A^/M^ ps 1, 
one needs F(l) = 14/3 — 2n and F(— 1) = 14/3 + 2n for a numerical estimate; with 
fi and g 2 as above, this term is 0.9% and —2.0% for D and D* , respectively. 

To put it differently, if the function F is expanded in powers of l/m c , then 
since the resulting terms have inverse powers of they cannot be confused 

with local counterterms. These terms play the same role as the non-analytic logs 
of the previous sections. 



= 2v, {^j 2 [F(A D /M n ) + ln(p 2 /M 2 )] +C(p)/m 2 j 

(6.19a) 

= K f 1 - Y (^7) 2 [F(-Ad/M v ) +Mv 2 /M 2 )] +C'(ri/m 2 c 



6.5. Trouble On The Horizon? 



Frequently the non-analytic corrections to relations that follow from the sym- 
metries are uncomfortably large. A case of much interest is the relation between 
the form factors f± and h for B — > K transitions, relevant to the short distance 
process b — > se + e~ , 

(K(p K ) | s^b | B{p B )) = f+ (p B + p K y + /_ (p B - p K Y , 
(K(p K ) I so-^b | B(p B )) = ih [(p B + PkY(pb - Pk) v - (pb + Pk) u {pb - PkY] , 

and the form factors for B — > 7rez/, 

(7f(p w ) i u^b i b( Pb )) = h (p B + p«r + f- (pb - p«y ■ 

In the combined large mass and chiral limits only one of these form factors is 
independent: 

m b h = f+ = -/_ = /+ = -/_ (6.21) 

In this limit, the ratio of rates for B — > Ke + e~ and B — > 7rez/ is simply given, 
in the standard model of electroweak interactions, by |Vts/Kib| 2 ) times a pertur- 
batively computable function of the top quark mass. If the relation (6.21) held 
to good accuracy one could thus measure a ratio of fundamental standard model 
parameters. 6 

The non-analytic, one-loop corrections to the relations in Eq. (6.21) have been 
computed. [56] The results are too lengthy to display here. Numerically, the violation 
to SU(S)v symmetry is found to be at the 40% level. 7 

The phenomenological lagrangian that we have been considering extensively 
neglects the effects of states with heavy-light quantum numbers other than the 
pseudoscalar - vector-meson multiplet. The splitting between multiplets is of the 
order of 400 MeV and is hardly negligible when one considers SU(S)v relations 
involving both 7r and K mesons. For example, consider the effect of the scalar - 

6 Another application of this relation has been discussed by I. Dunietz[55]. Assuming 
factorization in B — > ipX, ratios of CKM elements can be extracted from these two body 
hadronic decays. 

7 The large violation of SU(3)v symmetry affects as well the results of Dunietz (see 
previous footnote). 



pseudovector-meson multiplet. One can incorporate its effects into the phenomeno- 
logical lagrangian. To this end, assemble its components into a "superfield", akin 
to that in Eq. (6.1) for the pseudoscalar - vector multiplet: [57] 



S a (v) 



2 



K7m7 5 " B* 0a ] . 



(6.22) 



The phenomenological lagrangian has to be supplemented with a kinetic energy and 
mass for S, 



where A is the mass splitting for the excited S from the ground state H, and with 
coupling terms 



In terms of these one can now compute additional corrections to quantities such 
as /dJ/d in Eq. (6.10). Numerically the corrections are not small, [58] /d s //d = 
1 + 0.13/i 2 for M D * = 2300 MeV (or f D Jf D = 1 + O.O8/1 2 for M D * = 2400 MeV), 
assuming the strange mesons to be 100 MeV heavier. Similarly, corrections to the 
Isgur-Wise function can be computed, and are not negligible. [58] 
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Tr [S a {v)(iv ■ D ba - A8 ba )S b (v 
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